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6PSS  AND  MODELING  OF  COMPUTER  COMMUNICATION  NETWORKS 
1.  INTRODUCTION 

1.1  Objectives. 

In  order  to  determine  the  suitability  of  the  discrete  event 
simulation  language  GPSS  for  modeling  computer  network  structures  likely 
to  be  encountered  in  command,  control,  and  communication  (C3)  systems, 
several  example  computer  networks  were  simulated  using  this  language. 

This  report  presents  a  survey  of  GPSS  capabilities  and  peculiarities. 
Problems  encountered  in  translating  GPSS  programs  from  one  available 
version  of  GPSS  to  another  as  well  as  explanation  of  differences  in 
simulation  results  are  discussed.  Results  comparing  performance  of 
four  ring  network  structures  simulated  in  this  study  are  also  presented. 

1.2  Background. 

This  is  the  first  in  a  series  of  reports  on  the  progress  and 
results  of  AMSAA's  work  in  creating  new  and  using  existing  models  in 
analyzing  and  predicting  the  performance  of  computer  networks  and  their 
supporting  communications. 

1.2.1  Relation  to  Modeling.  The  study,  construction, 
and  validation  of  simulation  models  aids  the  work  of  the  group 
by  providing: 

•  data  to  support  conclusions  about  proposed  system  concepts, 

f  tools  for  evaluating  alternative  configurations  posed  by 
requirements  definition  studies  such  as  TOS  CASE  and  ASAS  FSD, 

•  the  means  to  examine  simultaneously  conputer  system  perform¬ 
ance,  network  configurations,  and  imperfect  communications, 

a  quantitative  estimates  of  the  effects  of  interoperability 
requirements  upon  the  performance  of  the  primary  mission 
of  a  system  and  upon  the  supporting  communications, 

•  data  to  augment  that  obtained  from  system  testing,  and 

•  estimates  of  the  most  difficult  combinations  of  system 
inputs  to  satisfy,  which  can  be  used  to  guide  test  planning 
toward  efficient  and  effective  discovery  of  system 
deficiencies. 

1.2.2  ^tiyation.  The  work  reported  here  was  motivated 
originally  by  an  effort  relating  to  TOS  CASE  in  which  varying  approaches 
were  proposed  by  contractors  for  modeling  and  simulating  the  combined 
computer  processing  and  communication  network  for  this  system.  Initially, 
a  model  of  conputer  communications  developed  for  the  Air  Force  called 
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SACDIN  was  proposed  for  use  but  subsequently  was  rejected  because  It 
could  not  easily  be  modified  to  handle  message  routing  in  other  than  a 
tree  connected  hierarchical  fashion.  The  contractor  then  proposed 
using  a  general  purpose  discrete  event  simulation  language  called  GPSS 
to  write  a  simulation  model  using  the  dialect  of  GPSS  implemented  on  a 
Control  Data  Corporation  6600  computer. 

To  prepare  AMSAA  personnel  for  analysis  of  the  validity  of  the 
anticipated  TOS  CASE  simulation  model,  a  study  of  GPSS  was  begun.  Because 
only  the  UNIVAC  dialect  called  GPSS  1100  v/as  available  to  AMSAA  personnel 
at  the  inception  of  this  study,  the  question  of  syntactically  and  seman¬ 
tically  correct  translation  of  simulation  programs  (models)  from  one  dialect 
of  GPSS  to  another  was  raised.  Much  of  the  work  reported  here  deals  with 
answering  this  question.  ^ 

1.2.3  Approach.  In  order  to  develop  expertise  in  GPSS  modeling 
of  computer  communication  networks  and  to  develop  confidence  in  comparison 
of  models  written  in  one  dialect  of  GPSS  with  those  written  in  a  different 
dialect,  the  team  decided  to  translate  known  computer  communication 
network  models  from  one  set  of  syntax  and  semantics  into  the  other. 

Models  of  computer  communications  networks  written  in  GPSS  were 
obtained  from  the  open  literature.  Only  those  having  a  relatively  simple 
structure  coupled  with  published  simulation  results  for  comparison 
purposes  were  considered  suitable  for  use  in  the  study  of  translation 
from  one  dialect  to  another.  Three  of  the  models  were  written  in  an 
enhanced  version  of  GPSS/360  for  an  IBM  360  series  computer.  The  fourth 
model  had  been  written  in  GPSS  1100  for  a  UNIVAC  1108  computer. 

Initially,  the  study  team  was  restricted  to  using  only  a  UNIVAC 
1108  conputer;  so  three  of  the  programs  were  translated  from  GPSS/360 
into  GPSS  lino,  and  several  differences  in  output  results  were  noted. 

Because  the  syntax  of  GPSS  1100  differs  from  that  of  GPSS/360  and  its 
later  dialect  called  GPSS/V,  the  correctness  of  the  syntactic  translation 
was  studied.  Careful  desk  checking  of  the  translation  by  at  least  three 
independent  programmers  revealed  no  discernable  errors,  leading  to  a 
check  of  possible  semantic  differences. 

Semantic  differences  are  those  due  to  the  way  in  which  the 
simulation  command  interpreter  is  actually  executed.  If  the  interpreter 
is  written  in  a  high  level  language  (e.g. ,  FORTRAN),  the  differences 
may  be  due  to  the  manner  in  which  the  various  subprograms  are  compiled; 
if  the  interpreter  is  written  in  assembly  language,  the  differences  may 
be  due  to  different  hardware  characteristics  such  as  word  length. 

Because  the  simulators  both  rely  on  pseudo  random  number  gener¬ 
ators  to  generate  the  stream  of  random  events  according  to  assumed  prob¬ 
ability  distributions,  it  was  first  necessary  to  account  for  possible  dif¬ 
ferences  generated  here.  Because  of  differences  in  word  length,  the 
largest  representable  integers  in  the  two  systems  are  different.  Hence, 
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the  two  pseudo  random  number  generators  are  Inherently  different.  Ini¬ 
tially,  it  was  postulated  that  either  one  or  both  sets  of  pseudo  random 
number  generators  may  be  exhibiting  nonrandom  behavior.  To  check  this 
hypothesis  some  tests  of  randomness  were  performed  on  the  generators,  and 
these  tests  are  documented  In  Appendix  A.  Even  If  the  generators  are 
sufficiently  random,  semantic  differences  In  the  way  In  which  the  generated 
numbers  are  subsequently  used  may  be  the  cause  of  differences  In  the  output 
results.  Deterministic  and  Identical  tables  of  numbers  supplied  to  both 
simulation  dialects  to  guide  event  generation  and  flow  In  the  models, 
coupled  with  detailed  traces  of  activity  In  the  models  were  then  considered 
appropriate  for  finding  differences.  This  approach  required  the  avail¬ 
ability  of  an  IBM  360  computer  system  or  equivalent.  Ultimately,  access 
to  an  IBM  360  computer  system  with  GPSS/360  was  obtained.  The  pseudo 
random  number  generators  In  both  the  GPSS/360  and  GPSS  1100  simulations 
were  disabled,  and  Identical  tables  of  pseudo  random  numbers  (generated  on 
a  CDC  7600)  were  appropriately  formatted  and  Inserted  Into  the  two  dif¬ 
ferent  dialect  simulation  models.  As  a  result,  certain  semantic  differences 
have  been  Identified  and  are  discussed  In  this  report. 

1.3  Organization. 

Chapter  2  of  this  report  discusses  briefly  the  concepts  of  dis¬ 
crete  event  simulation  and  presents  a  short  Introduction  to  the  GPSS  lan¬ 
guage  and  some  of  Its  capabilities  that  are  relevant  to  computer  communica¬ 
tions  network  simulations. 

Chapter  3  Introduces  the  ring  network  examples  simulated  In  this 
study  and  presents  results  of  those  simulations.  Lessons  learned  are  also 
discussed. 


The  appendices  Include  a  summary  of  pseudo  random  number  generator 
tests  and  their  results,  listings  of  GPSS  programs  for  the  computer  networks 
used  In  this  study,  and  a  glossary  of  acronyms  and  abbreviations. 

1.4.  Summary  of  Conclusions. 

Several  conclusions  were  reached.  Thqy  are; 

(1)  It  Is  possible  to  correctly  translate  simulation  programs 
from  one  dialect  of  GPSS  Into  another,  even  though  GPSS/360  and  GPSS  1100 
differ  In  both  syntax  and  semantics.  The  GPSS/360  syntax  uses  fixed  fields, 
and  GPSS  1100  has  a  column  free,  easier  to  use  format.  Semantic  differences 
are  due  to  Inherently  different  pseudo  random  number  generators,  the  docu¬ 
mented  use  of  differing  default  conditions,  and  undocumented  differences  In 
function  Interpolation.  Because  of  these  differences,  care  should  be 
taken  when  comparing  output  data  from  one  dialect  with  that  from  another. 

(2)  GPS$/360  on  the  APG  IBM  360/65  executes  considerably  faster 
than  does  GPSS  1100  on  the  ARRADCOM  UNIVAC  1108,  about  four  to  seven  times 
faster  for  the  examples  considered  here  and  for  other  test  cases  that  have 
been  run. 


(3)  Both  GPSS/360  and  RPSS  1100  have  attractive  features  for 
the  discrete  event  modeling  of  computer/communi cat ions  networks.  Messages 
are  easily  modeled  as  dynamic  entities  called  transactions.  Language 
features  are  provided  for  causing  message  arrivals  and  other  randomly 
occuring  events.  Equipment  entities,  such  as  transmitters,  receivers, 
and  message  queues  are  easily  modeled.  Automatic  collection  and  display 
of  statistics  on  system  performance  are  provided. 

(4)  Preliminary  analysis  of  end  to  end  message  transmission 
delay  data  from  simulations  of  four  ring  network  link  level  protocols 
indicates  that  at  low  system  loading  there  is  no  significant  (order  of 
magnitude)  difference  among  them.  The  systems  saturate  or  exhibit  expo- 
nentially  unbounded  end  to  end  delay  times  when  sufficiently  heavy  loads 
are  applied,  and  they  do  so  in  an  order  of  increasing  load  consistent 
with  previously  published  data. 

2.  SIMULATION  WITH  GPSS 

2.1  Discrete  Event  Simulation. 

System  simulation  using  models  having  state  variables  that  change 
state  only  at  discrete  instants  of  time,  with  time  progressing  in  dis¬ 
crete  increments,  is  called  discrete  event  simulation.  For  a  given 
interval  of  simulation  time,  points  of  event  occurrence  in  discrete  event 
simulation  are  both  finite  and  countable,  whereas  in  continuous  system 
simulation  the  time  of  event  occurrence  is  continuous.  Because  GPSS  is  a 
discrete  event  simulation  language,  any  system  being  modeled  in  GPSS 
must  be  representable  as  a  discrete  event  system.  Doing  so  requires 
what  may  appear  to  be  some  degree  of  oversimplification,  but  simpli¬ 
fication  is  acceptable  if  the  model  accurately  reflects  system  behavior 
without  necessarily  reproducing  exactly  and  completely  the  actual  system 
operation.  Since  there  are  many  different  simulation  languages  available 
to  the  user  the  features  that  distinguish  GPSS  will  be  examined. 

2.2  Features  of  GPSS. 

One  of  the  major  advantages  in  using  a  language  such  as  GPSS  to 
simulate  systems  is  the  convenience  afforded  by  the  language  [ll.  Instru¬ 
menting  a  simulation  model  to  collect  data  and  compute  statistics  reveal¬ 
ing  the  performance  of  system  components  of  interest  is  a  major  task  in 
constructing  a  system  model.  A  large  part  of  the  statistics  gathering 
is  intrinsic  to  GPSS;  hence  the  programmer  need  not  ordinarily  be  burdened 
with  this  time  consuming  task.  Along  with  its  automatic  data  collection, 
GPSS  allows  the  modeling  of  many  of  the  significant  characteristics  of 
"real  <yorld"  systems  with  much  ease.  The  characteristics  that  are  easily 
represented  include  dynamic  entities,  equipment  entities,  operational 
entities,  data  entities,  and  randomness  considerations.  Also  subliminally 
used  are  the  simulation  clock  and  the  event  scheduling  algorithm.  A 
brief  description  of  each  of  these  factors  follows. 
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2.2.1  Dynamic  Entities.  The  (tynamlc  entities,  called  transac> 
tions  in  GPSS,  are  used  to  represent  a  flow  of  some  sort  through  the 
system.  The  transactions  which  "flow"  through  the  model  may  either 
cause  an  activity  or  be  the  recipient  of  an  activity.  In  other  words, 
the  transaction  may  itself  cause  the  state  of  the  model  to  change,  or  it 
may  have  any  of  its  associated  parameters  changed.  The  altering  of  a 
parameter  value  of  a  transaction  may  in  turn  be  used  at  another  place 
(or  time)  in  the  model  to  effect  changes  to  the  state  of  the  model. 

2.2.2  Equipment  Entities.  Equipment  entities  are  used  in 
modeling  components  tnat  have  a  specific  action  associated  with  them. 
Equipment  entities  include  storages,  facilities,  and  logic  switches. 
Storages  are  used  to  represent  entities  that  may  have  their  activity 
dictated  by  one  or  more  transactions,  whereas  facilities  are  used  to 
represent  entities  that  may  have  their  activity  dictated  by  only  one 
transaction  at  a  time.  A  logic  switch  is  used  as  a  binary  state  indicator, 
such  as  locked  or  unlocked,  available  or  unavailable,  and  open  or  closed. 

2.2.3  Operational  Entities.  The  operational  entities  are 
used  to  perform  a  variety  of  functions.  They  provide  for  representation 
of  system  relationships,  model  activity  control,  and  the  basic  structure 
of  the  model  to  name  only  a  few.  In  SPSS  the  operational  entities  are 
blocks,  queues,  user  chains,  groups,  and  save  values.  Blocks  are  the 
basic  unit  of  the  model  structure.  Queues  are  generally  used  to  monitor 
delays  encountered  by  transactions  at  specific  points  in  the  model. 

User  chains  are  used  to  alter  the  normal  "flows"  of  transactions  in  a 
user  defined  manner.  Transaction  "flov/"  may  be  controlled  on  the  basis 
of  group  membership,  where  group  membership  indicates  a  certain  relation¬ 
ship  existing  between  transactions  in  the  group,  Savevalues  are  used  to 
store  information  at  certain  locations  in  the  model. 

2.2.4  Data  Entities.  Data  entities  are  used  to  input  data 
to  and  output  data  from  the  model  as  well  as  to  represent  certain  data 
relationshios.  The  data  entities  available  to  the  RPSS  user  include 
functions,  variables,  and  tables.  Functions  are  a  means  of  entering 
distributions  of  various  types  to  the  model.  The  distributions  may 
represent  real  system  data  or  they  may  merely  specify  standard  distribution 
forms  that  may  be  necessary  to  the  simulation.  Variables  are  used  to 
represent  system  data  relationships.  Tables  are  included  as  a  means  of 
extracting  data  from  the  model. 

2.2.5  Pseudo  Random  Number  Generators.  In  addition  to  the 
various  entities  that  can  be  modeled,  the  GP$S  programmer  has  a  number 
of  pseudo  random  number  generators  available  to  him  to  aid  in  the  simula¬ 
tion  of  randomly  occurring  events.  The  pseudo  random  number  generators 
are  actually  deterministic,  of  course,  but  this  offers  one  distinct 
advantage--reproducibility  of  simulation  runs  for  program  debugging 
purposes  using  the  same  sequence  of  numbers  from  run  to  run. 

2.2.6  Simulation  Clock  and  Event  Scheduling  Algorithm.  The 

s i mu 1  a t i on  clock  and  the  event  scheduling  algorithm  are  related  concept s . 
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The  GPSS  simulation  clock  does  not  advance  time  In  fixed  unit  Increments. 
Instead  the  simulation  clock  Is  advanced  only  when  the  next  event  Is 
scheduled,  and  Is  advanced  to  that  next  scheduled  time  directly.  Event 
scheduling  Is  effected  by  scanning  one  of  several  "event  chains,"  or 
ordered  lists  of  transactions.  After  the  appropriate  chain  Is  scanned, 
processing  of  transactions  occurs  on  the  basis  of  scheduled  departure 
time,  currently  assigned  priority,  and  time  resident  on  the  chain. 

After  all  events  that  can  take  place  at  the  current  simulation  clock 
time  have  occurred,  the  simulation  clock  Is  advanced  to  the  next  scheduled 
event  occurrence  determined  by  a  scan  of  the  future  events  chain.  Simula¬ 
tion  continues  In  this  fashion  until  an  event  occurs  that  terminates  the 
simulation  at  some  desired  point. 

2.3  Comparison  of  GPSS/360  and  GPSS  1100. 

The  two  dialects  of  GPSS  available  to  the  study  team  were  IBM's 
6PSS/360  [2]  and  UNIVAC's  GPSS  1100  [31.  The  IBM  version  of  GPSS  executes 
on  the  APG  IBM  360/65  computer  system,  and  the  UNI VAC  version  executes  on 
the  ARRAOCOM  UNIVAC  1108  computer  system.  These  two  versions  of  GPSS  are 
distinct  Implementations  of  the  same  discrete  event  simulation  concept, 
but  there  are  a  number  of  differences  between  them  as  discussed  below. 

2.3.1  Syntax.  Both  versions  of  GPSS  have  the  same  basic 
block  structure,  but  syntax  varies  considerably  between  the  two.  UNIVAC 
GPSS  1100  uses  a  relatively  free  form  Input  format  In  Its  statement 
specification  language.  Similar  to  the  UNIVAC  Assembler  Input  statement 
formats,  various  fields  appearing  In  the  line  Image  of  a  GPSS  1100  source 
statement  are  not  column  dependent,  are  simply  separated  by  one  or  more 
blank  spaces,  and  in  some  cases  are  not  required  to  appear  In  a  specified 
order.  In  IBM  GPSS/360  the  fields  of  a  source  statement  must  appear  In 
specific  column  locations  In  the  line  Image.  For  example,  the  location 
field  used  to  Identify  a  specific  statement  for  later  symbolic  reference 
must  begin  1n  column  two  and  not  extend  past  column  six.  This  places  a 
five  character  limitation  on  statement  names  or  Identifiers.  Identifiers 
In  GPSS  1100  can  be  more  than  five  characters  In  length,  resulting  In 
the  ability  to  use  more  descriptive  location  names. 

2.3.2  Function  Definition.  Another  difference  between  the 
languages  Is  in  the  area  of  user  defined  versus  simulator  supplied  func¬ 
tions.  Both  simulators  provide  several  callable  pseudo  random  number 
generators  with  which  simulator  supplied  uniform  distribution  functions 
are  generated  and  for  v/h1cb  the  user  need  only  supply  the  end  points.  In 
order  to  specify  an  exponential  probability  distribution  function  or  a 
Gaussian  distribution  function  In  IBM  GPSS/360,  the  user  must  supply  a 
finite  set  of  x  and  y  coordinates  that,  coupled  with  simulator  supplied 
linear  Interpolation,  approximate  the  desired  distribution  function. 
Oepending  on  desired  accuracy,  approximations  of  24  to  60  or  more  points 
are  typically  used.  The  UNIVAC  GPSS  1100  simulator  supplies  uniform, 
exponential,  and  Gaussian  distribution  functions  as  built-in  components 
of  the  language.  To  Invoke  the  exponential  or  Gaussian  distribution, 
the  GPSS  1100  user  need  only  reference  them  with  appropriate  parameters 
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(a  mean  value  for  the  exponential  function  and  a  mean  and  variance  for 
the  Gaussian  function).  As  with  the  IBH  GPSS  Implementation,  the  user 
can  define  any  other  desired  functions  by  specifying  appropriate  sets  of 
points. 


2.3.3  Memory  Allocation.  Because  IBM  GPSS/360  allows  the 
programmer  to  speciry  eitner  nairword  or  fullword  values  for  parameters 
and  savevalues,  the  programmer  can  save  some  memory  space  for  allocation 
to  other  purposes.  This  represents  an  advantage  over  the  GPSS  1100 
version.  The  assignment  of  halfword  parameters  and  savevalues  normally 
might  be  used  only  minimally  by  most  modelers.  A  second  and  fairly 
small  benefit  Is  that  smaller  models  run  faster.  Perhaps  there  are  only 
a  few  Instances  where  a  decreased  run  time  may  be  noticeable,  but  In 
these  few  Instances  It  may  be  a  large  advantage.  The  feature  of  variable 
word  size  for  parameters  and  savevalues  gives  GPSS/360  greater  flexibility 
than  GPSS  lion. 


2.3.4  Function  Interpretation  (Interpolation).  The  two  versions 
of  GPSS  differ  slightly  in  the  way  that  they  perform  interpolation  In' 
user  defined  functions.  For  example, a  continuous  function  may  be  defined 
with  x-coordinate  values  of  0  and  1000  and  corresponding  y-coordinate 
values  of  1  and  6,  respectively.  This  defines  a  straight  line  segment 
between  the  points  (0,1)  and  (1000,6).  Now,  given  that  the  x  value  Is 
to  be  determined  by  some  random  number  generator  with  values  ranging 
from  0  to  999,  and  that  both  Interpreters  operate  by  truncation  rather 
than  rounding,  the  functions  can  then  yield  results  of  1,2, 3, 4,  or  5 
with  equal  likelihood.  Since  the  representation  of  s1ngle>prec1s1on 
floating  point  numbers  In  IBM  360  computers  uses  a  32-b1t  hexadecimal ly 
normalized  format  and  In  UMIVAC  1100  computers  a  36-b1t  binary  normalized 
format,  the  representation  of  certain  fractions  Is  not  exact. 

The  expression  of  certain  nur4»ers  was  found  to  be  a  problem  In 
the  above  example.  It  was  found  that  for  an  x  value  of  200,  the  IBM 
simulator  returned  a  y  value  of  2~the  result  that  one  would  expect. 

The  UHIVAC  simulator,  however,  returns  a  value  of  1  for  the  same  Input 
X  value.  Further  Investigation  found  that  both  the  IBM  and  UNIVAC 
versions  returned  the  correct  y  value  of  2  for  an  x  value  of  201,  and 
the  correct  y  value  of  1  for  an  x  value  of  199.  The  problem  again  arose 
In  the  evaluation  of  x  coordinates  of  400,  6no,  and  800. 

One  reason  for  the  discrepancy  may  be  attributed  to  the  order 
In  which  arithmetic  operations  are  carried  out  In  the  Interpolation 
process.  Since  truncation  Is  used,  the  order  of  operations  is  Important. 

For  example,  letting  (xi,yi)  and  (X2,y2)  be  the  endpoints  of  a  continuous 
stralghtllne  function  In  which  Intermediate  Interpolated  values  are 
desired,  the  Interpolated  value  y  Is  given  by: 

y  ■  r(y2-yi)/(x2-Xj)  ]  -x  +  r(x2yi  -  Xjy2)/  (x2-xi)] 


*  mx  +  b 


where 


•n  ■  C(y2-yi)  /  (X2-Xi)1  .  and 

b  ■  yj  if  xj  «  0. 

In  the  case  considered  here,  b  ■  yj  . 

Two  of  the  possible  combinations  for  ordering  operations  in  the 
computation  of  y  are: 

Approach  1: 

Step  1:  set  m  :»  [(y2-yi)/(x2-xi )1 
Step  2;  set  z  :»  m  •  x 
Step  3:  set  y  :*  2  +  b 

Step  4:  set  y  integer  [y]  ,  i.e.  truncate  fraction. 

Approach  2: 

Step  1:  Set  z  :»  (y2-yi)*x 
Step  2:  Set  w  :»  z/(x2-xi) 

Step  3:  Set  y  :»  w  +  b 
Step  4:  Set  y  :»  integer  Cy1  . 

In  certain  instances  such  as  (xi,yi)  ■  (0,1)  and  (x2,y2)  • 
(1000,6)  and  x  =  200,  Step  3  of  Approach  1  produces  1.9999999926iQfor 
the  UNIVAC  single  precision  floating  point  format  and  1, 99999904b jiq 
for  the  IBM  single  precision  floating  point  format.  If  the  order  of 
operations  in  both  IBM  and  UNIVAC  GPSS  implementations  corresponds  to 
Approach  1  (and  at  least  IBM  GPSS/360  documentation  [221  pp.  75  &  205 
seems  to  so  indicate),  then  the  y  value  returned  in  both  systems  (after 
truncation)  would  be  unity.  Using  Approach  2  with  the  same  data  items 
as  above,  the  result  is  the  integer  value  y«2  for  both  the  UNIVAC  and 
and  IBM  interpolation  schemes.  Empirical  results  using  the  above  data 
items  in  both  GPSS  implementations  produces  interpolated  integer  values 
of  y  “  1  for  the  UNIVAC  inplementation  and  y  ■  2  for  the  IBM  implement- 
ation,  indicating  that  perhaps  the  available  IBM  GPSS  documentation 
does  not  accurately  reflect  the  actual  ordering  of  operations,  or  that 
the  documentation  available  to  the  study  team  does  not  include  all 
possible  change  notices.  The  UNIVAC  inplementation  would  appear  from 
this  single  sample  to  accurately  follow  the  operation  ordering  stated 
in  the  IBM  documentation.  In  any  event  a  likely  cause  of  observed 
differences  in  GPSS  function  interpolation  between  the  two  implementa¬ 
tions  is  due  to  different  (nonequivalent)  orderings  of  finite  precision 
floating  point  arithmetic  operations. 


Determining  the  exact  cause  of  the  differences  would  require 
laborious  and  time  consuming  detailed  examination  of  the  assembly  level 
machine  code  for  the  two  GPSS  Implementations,  which  Is  beyond  the  scope 
of  this  stu<1y.  The  most  Inportant  fact  has  been  ascertained:  namely, 
exact  and  correct  syntactic  translations  of  GPSS  programs  between  IBM 
GPSS/360  and  UNIVAC  GPSS  1100  can  produce  differing  output  values  that 
are  caused  by  semantic  differences  In  the  Implementations  of  Interpolation. 

Differences.  Miscellaneous  differences 
between  GPSs/36d  and  GPsS  lloO  Include  the  sinulatlon  clock  starting 
time  and  the  calculation  of  standard  deviations  In  the  standard  statisti¬ 
cal  output.  The  IRM  version  of  GPSS  starts  Its  simulation  clock  at  time 
one,  while  the  UNIVAC  version  starts  Its  simulation  clock  at  time  zero. 

This  is  a  minor  difference,  hut  one  whose  effect  can  be  seen  when  a 
model's  transaction  routing  is  a  function  of  absolute  simulation  clock 
time.  The  UNIVAC  clock  can  be  aligned  with  the  IBM  clock  by  specifying 
that  no  transaction  enter  the  model  before  time  one.  Differences  In 
calculated  standard  deviations,  though  tmall,  were  observed  when  start 
time,  and  the  generation  and  movement  of  all  transactions  were  forced  to 
be  identical  In  deterministic  models.  The  reason  for  these  standard 
deviation  differences  Is  not  apparently  due  to  one  version  producing 
best  estimates  of  standard  deviation  and  the  other  not  doing  so,  and  the 
exact  reasons  for  these  modest  differences  are  not  yet  understood. 

2.3.6  Random  Processes.  One  point  to  be  considered  when 
running  stochastic  simulations  ^s  whether  processes  to  be  modeled  as 
random  can  be  modeled  acceptably.  Each  of  the  two  versions  of  GPSS 
offers  pseudo  random  number  generators  to  aid  the  modeling  of  stochastic 
processes.  IBM  GPSS/360  offers  one  such  generator  replicated  eight 
times.  Hence,  a  user  can  Implement  up  to  eight  distinct  sequences  of 
random  numbers.  The  sequences  will  be  Identical  Initially  unless  the 
user  Inputs  a  seed  different  from  the  default  value  to  one  or  more  of 
the  generators.  UNIVAC  GPSS  1100  offers  ten  distinct  pseudo  random 
number  generators.  The  generators  are  of  the  same  type  (either  linear 
or  mixed  linear  congruentlal)  but  use  different  seeds  and  multipliers. 
Statistical  properties  of  pseudo  random  number  generators  for  both  GPSS 
versions  were  studied  to  determine  whether  the  generators  are  random 
enough,  and  details  of  that  study  are  presented  In  Appendix  A.  In  summary 
the  pseudo  random  number  generators  are  generally  random  enough  for  use 

In  the  ring  network  simulations  discussed  In  the  next  chapter. 

2.3.7  Run  Time.  One  last  consideration  of  the  differences 
between  IBM  GPSS/36fl  and  UNIVAC  GPSS  1100  Is  simulation  execution  time 
(or  run  time)  and  Its  corresponding  cost.  The  CPU  time  for  four  ring 
network  models  using  the  IBM  GPSS  simulator  was  from  one  fourth  to  one 
tenth  of  that  required  to  execute  the  same  models  using  the  UNIVAC  GPSS 
simulator.  For  example,  GPSS  simulation  of  the  DLCN  model  described  In 
Chapter  3  required  4  min.  16  sec.  of  CPU  time  for  the  IRM  version  and 
30  riiln.  33  sec.  of  CPU  time  for  the  UNIVAC  version  of  the  model  using 
Identical  system  parameters.  For  this  example  the  UNIVAC  version  runs 
about  seven  times  slower  than  the  IBM  version.  There  Is  apparently  a 
significant  speed  (and  hence,  cost)  advantage  in  running  GPSS/360  models 
over  the  GPSS  1100  models. 
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Turnaround  time,  measured  using  wall  clock  time,  was  also  generally 
better  on  the  APG  IBM  360/65  than  on  the  ARRADCOH  UNI VAC  1108  when  running 
corresponding  GPSS  simulations  for  the  four  ring  networks  considered  In 
Chapter  Three.  Mall  clock  time  Includes  a  measure  of  system  congestion, 
and  to  the  programmer  fast  turnaround  Is  usually  of  Interest.  Sample 
simulations  run  as  the  only  batch  job  on  the  system  at  times  when  time 
sharing  demand  service  was  cut  off  Indicate  similar  ratios  of  wall  clock 
time.  Sobe1[7]  was  plagued  by  extraordinarily  long  run  times  under  similar 
system  loading  conditions  on  a  UNIVAC  1100/42  system.  Simulation  runs 
that  finished  normally  on  the  APG  IBM  360/65  In  an  hour  of  wall  dock  time 
terminated  abnormally  on  the  much  faster  UNIVAC  1100/42  system  In  approxi¬ 
mately  four  hours  of  wall  clock  time  on  an  essentially  empty  system,  where 
abnonnal  termination  was  caused  by  the  need  to  exceed  the  programmer  speci¬ 
fied  run  time  limit.  Although  the  UNIVAC  1108  Is  a  faster  computer  than 
the  IBM  360/65  according  to  Schrlber  [1]  the  UNIVAC  GPSS  1100  simulator 
appears  to  have  a  far  less  efficient  Implementation  than  does  the  IBM 
GPSS/360  simulator.  Models  executed  from  four  to  seven  times  faster  In 
the  IBM  version.  In  addition,  comparison  of  wall  clock  times  for  the  four 
ring  network  simulations  revealed  that  the  IBM  360/65  system  gives  from 
two  to  three  times  better  turnaround  than  does  the  UNIVAC  1108  system. 

This  may  not  be  true  In  all  modeling  situations,  but  for  the  rather  simple 
ring  network  structures  studied  IBM  GPSS/360  is  more  efficient  than  UNIVAC 
GPSS  1100.  This  conclusion  Is,  of  course,  system  configuration  and  site 
dependent. 


2.4  Suitability. 

2.4.1  Ease  of  Model  Implementation.  The  first  factor  In  deter¬ 
mining  the  suitability  of  GPSS  for  modeling  conputer  communications  networks 
Is  the  ease  of  model  Implementation.  Each  block  In  the  structure  of  a 

GPSS  model  may  represent  a  separate  action  block  In  a  flowchart  of  the 
system  being  modeled.  For  Instance,  the  process  of  capturing  a  facility 
for  some  length  of  time  and  then  relinquishing  control  of  the  facility 
requires  three  GPSS  blocks:  one  to  seize  the  facility,  one  to  advance  the 
clock,  and  one  to  return  the  facility  to  Its  previous  state.  This  is 
considerably  simpler  to  specify  In  GPSS  than  it  might  be  In  many  other 
programming  languages  In  which  It  may  be  necessary  to  write  one  routine  to 
implement  each  of  the  three  GPSS  blocks.  The  event  processing  routines 
are  Intrinsic  to  the  GPSS  language,  so  the  user  need  not  be  bothered  by 
the  possibly  unpleasant  task  of  describing  each  action  In  detail. 

2.4.2  Understandablllty.  Another  factor  in  the  ease  of  model 
Implementation  In  GPSS  is  this  language's  choice  of  block  names  which  aids 
understandablllty.  The  process  of  obtaining  control  of  some  facility  Is 
written  as  SEIZE  "facility"  In  GPSS.  The  SEIZE  block  Is  then  a  model 
statement  that  can  be  readily  understood  by  managers  as  well  as  program¬ 
mers.  The  majority  of  blocks  In  GPSS  are  named  in  such  a  way  that  the 
block  name  describes  the  block  function. 
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2.4.3  Standard  Statistical  Output.  Another  advantaqe  In 
building  models  with  01*55  is  the  standard  statistics  gathering  Intrinsic 
to  and  aided  hy  the  language.  Statistics  such  as  gueiie  times,  storage 
contents,  and  facility  utilizations  are  all  collected  automatically  by 
the  OPSS  simulator.  These  Items,  along  with  a  large  number  of  other 
useful  statistics,  are  printed  In  a  standard  statistical  package  In  the 
output  report  of  the  program.  Additional  Information  concerning  the 
model  run  can  be  obtained  by  the  Inclusion  of  user  defined  tables  In  the 
output  report. 

2.4.4  Optional  Output.  As  optional  output,  TRACE  and  PRINT 
blocks  are  available  to  aid  In  the  debugging  of  a  GPSS  program.  After 

all  known  bugs  have  been  removed  from  the  simulation  model,  the  programmer 
may  specify  optional  output  formats  and  histograms  as  well  to  make  the 
output  understandable  to  the  np*?specl8l1st. 

2.4.5  Level  Pgtall.  An  additional  consideration  In  assess- 
Inn  the  appropriateness  of  RPSS  for  computer  communication  network  models 
is  the  level  of  detail  permitted  In  the  models.  If  the  modeling  ob.lective 
is  to  develop  an  exact  detailed  replica  of  the  real  world  system,  then 

It  is  doubtful  that  GPSS  would  be  a  suitable  language.  If,  however,  the 
objective  of  the  model  Is  to  gain  general  Insights  Into  how  a  system 
will  perform  under  various  circumstances,  then  GPSS  could  he  a  suitable 
language.  Because  there  are  memory  space  limitations  on  the  size  of  the 
GPSS  program,  some  si ppHTI cations  must  be  made  as  a  trade-off.  In 
deciding  whether  to  model  In  GPSS,  the  analyst  must  determine  whether 
the  amount  of  simplification  regulred  is  acceptable.  Language  features 
pemit  the  programmer  to  command  reallocation  of  the  available  data 
storage  space  among  the  competing  entitles  Invoked  iy  block  specifications. 
However,  large  models  (I.e.,  those  with  large  numbers  of  blocks  or  large 
numbers  of  sirvltaneously  active  transactions)  can  easily  exceed  the 
available  storage  on  the  machine  executing  the  GPSS  simulator.  In  such 
cases  the  programmer  may  he  forced  to  redtice  the  level  of  detail  simulated 
In  order  to  get  his  model  to  run  at  all  in  the  existing  har«Vare/software 
environment. 

Similar  decisions  and  limitations  are  faced  by  analysts  and 
programmers  In  every  language  chosen  for  performino  simulation.  In  some 
lanouanes  the  ability  to  call  operating  system  service  rotitines  or  other 
library  routines  may  be  more  easily  performed  than  In  GPSS.  Resolving 
problems  at  acceptable  cost  in  tim**  and  effort  is  the  key  Issue  and 
must  be  traded  against  ease  of  simulation  model  Implementation  directly 
available  from  language  features  and  level  of  detail  required. 

3.  COMPUTER  COMMUNICATION  NETWORK  MODELS 

3.1  Network  Concepts  and  Terminology. 


Computer  communication  networks  are  essential  components  of 
military  c3  systems.  Computer  communication  networks  permit  users  to 
access  resources  such  as  har(bfare  units,  software  packages  and  data  files 
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In  a  remote  computer  system.  One  can  view  the  structure  of  a  computer 
communication  network  as  heino  partitioned  Into  two  parts,  a  communi¬ 
cation  network  (sometimes  called  the  communication  subnetwork)  and  a 
user  resource  network (41. 

3.1.1  Communication  Network.  The  communication  network  comprises 
the  switchinq  computers  (or  nodes)  and  the  communication  channels.  Its 
function  Is  to  deliver  messaqes  from  one  node  to  another. 

3.1.2  User  Resource  Network.  The  collection  of  terminals  and 
computinq  resources  comprises  the  user  resource  network.  These  resources 
are  connected  to  switchinq  nodes  and  communicate  with  each  other  by  way 
of  the  communication  network. 

3.1.3  Hosts.  Protocols,  and  Network  Function.  The  computer 
systems  In  the  user  resource  network  are  called  hosts,  ar1  a  set  of 
protocols  Is  Implemented  In  the  operatinq  system  of  each  host.  These 
protocols  are  procedures  to  Initiate,  maintain  and  terminate  software 
communications  via  the  nodes  of  the  communication  network.  A  host 
computer  may  accent  jobs  (such  as  requests  for  processing,  data  base 
queries  or  updates,  etc.)  from  local  or  remote  users.  Remote  jobs  are 
received  as  messaqes  from  the  communication  network,  and  require  extra 
processing  time  for  orotocol  handling.  When  orocessinq  of  the  remote 
task  Is  copmlete,  the  results  are  repackaged  as  a  message  (or  a  set  of 
related  messages)  and  are  returned  to  the  remote  users  via  the  communi¬ 
cation  network. 

3.1.4  Message  Switching.  The  basic  technique  hv  which  messages 
are  dellvefW  ft*W  Sobrce  node  to  destination  node  In  a  communication 
network  Is  called  message  switching.  In  this  technique  a  message  entering 
the  network  Is  first  passeo  to  its  origin  node  where  It  may  be  stored 
while  It  waits  for  route  selection  according  to  some  routing  algorithm 
and  where  It  mav  queue  for  Its  outbound  communication  channel.  When  the 
channel  becomes  free,  the  message  Is  transmitted  to  the  next  node  along 
Its  route  to  the  destination,  and  the  above  process  Is  repeated  until 

the  message  Is  delivered  to  Its  destination  node. 

3.1.5  Packet  Switching.  '  A  modification  of  message  switching 
Is  a  technique  called  packet  switching  wherein  each  message  is  decomposed 
Into  maximum  length  disjoint  subsets  called  packets.  Each  packet  Is 
Identified  for  later  message  reassembly,  and  each  can  be  routed  Independ¬ 
ently  through  the  communication  network. 

3.1.6  Performance  Measures.  The  total  elapsed  time  from  the 
arrival  of  a  message  at  its  source  node  to  the  successful  delivery  of 
this  message  to  Its  destination  Is  called  end-to-end  delay  and  Is  an 
Important  performance  measure  of  both  message  and  packet  switched  networks. 
Factors  Influencing  this  performance  measure  Include  assumed  (or  actual) 
message  arrival  and  message  length  statistics,  routing  algorithms,  channel 
service  and  error  rates,  resource  contention  and  assigned  priority  classes, 
and  queueing  and  buffering  delays  enroute.  In  order  to  minimize  end-to-end 


delay,  designers  need  tools  with  which  to  predict  Its  mean,  variance, 
and  distribution  subject  to  sets  of  Input  parameters.  Other  performance 
measures  and  the  effects  of  design  parameters  nust  also  be  analyzed  In 
order  to  determine  quantities  such  as  optimal  finite  buffer  size,  channel 
utilization,  and  system  throughput  (I.e.,  messages/unit  time). 

3.2  Network  Modeling. 


Queueing  network  models  have  been  used  extensively  In  the  per¬ 
formance  analysis  of  message  switched  (or  packet-switched)  communication 
networks. 


3.2.1  Analytic  Models.  Closed  form  analytic  models,  when 
available,  are  advantageous  in  mat  they  can  lead  to  low  computational 
cost  predictions.  Exact  analytic  analyses  are  restricted  to  certain 
classes  of  simplified  models  [S],  and  results  for  general  models  with 
more  complex  features,  such  as  adaptive  routing  algorithms  and  finite 
buffer  space,  are  not  yet  available. 

3.2.2  Simulation  Models.  Discrete  event  simulation  models 
have  been  used  botn  to  verity  tne  adequacy  of  simplified  analytic  models 
and  to  provide  performance  analyses  In  cases  as  yet  too  complex  for 
adeguate  analytic  models.  The  generality  of  simulation  models  Is  paid 
for  In  higher  computational  costs  and  generally  greater  computer  execution 
times  than  may  be  required  for  evaluating  analytic  models.  If  partial 
analytic  results  are  available,  mixed  analytic  and  simulation  models 

help  to  reduce  simulation  costs.  In  many  cases  the  system  description 
parameters  such  as  non-Poisson  arrival  statistics  and  state  transition 
probabilities  are  either  not  available  or  not  directly  useable  In  the 
analytic  models;  whereas  enough  Information  may  be  available  to  Implement 
a  discrete  event  simulation  whose  Input  Is  a  list  of  measured  arrival 
events  from  some  actual  systems. 

3.3  Network  Topologies. 

Figure  3.1  shows  three  basic  topologies  commonly  found  In  com¬ 
puter  comrxinicatlon  networks;  the  mesh,  the  tree,  and  the  ring;  vari¬ 
ations  of  these  also  commonly  occur.  Internetwork  configurations  wherein 
nodes  in  one  topological  network  structure  act  as  gateways  to  other  (or 
even  the  same)  topological  network  structures  are  also  frequently  encountered. 

3.3.1  Mesh.  A  mesh  connection  of  nodes  Is  characterized  by  a 
connectivity  generally  greater  than  or  equal  to  two  at  each  node  so  that 
at  least  a  subset  of  nodes  can  select  alternate  routing  paths  between 
source-destination  pairs. 

3.3.2  Tree.  A  tree  connection  Is  characterized  by  a  hierarch¬ 
ical  structure  In  wnich  the  message  path  between  two  nodes  at  the  same 
level  In  the  tree  must  pass  through  a  common  ancestor  node  at  a  higher 
level  In  the  tree. 
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Figure  3.1:  Some  Computer  Commuuicat i jn  N'etwork  Topologies 
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3,3.3  Ring.  The  ring  Is  characterized  by  a  node  connectivity 
of  exactly  two  and  a  unidirectional  transfer  of  information  around  the 
comnunication  links.  A  message  going  from  a  given  node  to  its  predecessor 
node  in  the  ordering  inplied  by  the  direction  of  information  transfer  on 
the  ring  must  pass  through  all  the  nodes  on  the  ring  to  reach  its  destina¬ 
tion. 


3.3.4  Variations.  Topological  variations  in  mesh  connections 
range  from  minimal  to  maximal  connectivity,  and  to  structures  resembling 
tree  structures  with  cross  connections  between  a  subset  of  nodes  in 
different  branches  but  at  the  same  level  in  the  tree.  The  principal 
topological  variation  for  ring  (or  loop)  networks  conprises  two  or  more 
rings  (usually  passing  messages  in  opposite  directions)  for  greater 
reliability  and  increased  throughput. 

3.4  Ring  Network  Structures  Considered. 

Because  the  routing  structure  of  rings  and  loops  is  determin¬ 
istic  and  sinple,  and  because  GPSS  models  of  both  message  switched  and 
packet  switched  ring  networks  are  readily  available  in  the  literature 
[6,7],  this  network  topology  was  chosen  for  further  investigation  in  the 
simulation  study  of  computer  communication  networks  presented  here. 
Validation  of  the  simulation  models  and  comparisons  with  prior  work  of 
others  are  possible  for  this  topology  because  earlier  simulation  results 
are  available  in  the  open  literature  C8,9],  and  this  provides  greater 
dociimentation  and  insights  than  are  usually  available  for  more  complex 
topological  structures. 

A  loop  network  is  sometimes  distinguished  from  a  ring  network 
according  to  whether  the  communication  access  control  protocol  is 
centralized  (loop)  or  distributed  (ring).  Some  authors  refer  to  loops 
and  rings  interchangeably,  including  those  who  have  designed  loop  networks 
with  distributed  control  mechanisms  [8,9,10,11,12]. 

Four  basic  types  of  single  loop  networks  have  been  proposed  in 
the  literature,  namely,  the  Newhall,  Pierce,  DLCN,  and  Playthrough  struc¬ 
tures.  These  loop  networks  are  distinguished  by  their  transmission 
control  and  link  access  mechanisms, 

3.4.1  Tlie  earliest  loop  structure  was  proposed  by 

Farmer  and  Newhall  [13.1,  and  is  commonly  referred  to  as  a  Newhall  loop. 

In  this  structure  a  single  control  token  is  passed  from  one  loop  interface 
to  the  next  until  it  reaches  a  node  with  a  message  to  transmit.  That 
node  temporarily  seizes  the  control  token  and  starts  transmitting  its 
(variable  length)  message  to  an  addressed  destination  node.  Intervening 
nodes  pass  this  message  to  the  destination  node  which,  according  to 
varying  inplementations,  either  copies  the  message  into  its  arrivals 
buffer  or  removes  the  message  from  the  loop.  For  error  checking  purposes 
the  message  sometimes  is  permitted  to  circulate  to  the  receiver  portion 
of  the  source  node,  which  then  performs  a  consistency  check  and  removes 
the  message  from  the  loop.  Also,  depending  upon  implementation,  the 
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source  node  currently  In  possession  of  the  control  token  may  transmit 
one  or  more  variable  length  messages  before  relinquishing  control  of  the 
loop  by  passing  the  control  token  to  the  next  node  in  sequence.  Only  one 
source  node  may  transmit  at  a  given  time,  and  all  other  potential  source 
nodes  ntjst  wait  to  transmit  queued  messages  until  they  receive  the  control 
token.  Several  experimental  and  commercial  loop  communications  systems 
for  interconnecting  computers  and  conponents  have  been  based  on  minor 
variations  of  this  link  level  protocol  structure  (e.g.,  [14,  15]). 

3.4.2  Pierce.  The  Pierce  loop  [16,17,18,191  divides  conminica- 
tion  space  on  the  loop  into  an  integral  number  of  fixed  size  slots, 
called  packet  frames,  into  which  data  packets  can  be  placed.  To  send  a 
message,  a  node  segments  the  message  into  fixed  length  packets,  appends 
necessary  overhead  information  to  identify  both  the  packet's  number  and 
the  message  to  which  it  belongs,  places  each  packet  into  the  next  available 
empty  slot  passing  the  node,  and  marks  the  slot  as  full.  As  this  full 
message  packet  proceeds  toward  its  destination,  the  other  nodes  along 
the  route  examine  the  header  information  in  each  packet  frame  to  ascertain 
which  of  them  is  the  addressed  destination.  The  destination  node,  having 
recognized  its  address,  copies  the  data  being  received  and  either  fills 
the  slot  with  new  outbound  information  or  passes  this  now  empty  slot  to 
the  next  node.  Incorporated  into  the  loop  is  a  single  special  control 
node  that  maintains  time  slot  synchronization  for  the  loop  and  prevents 
buildup  of  undeliverable  packets.  The  header  of  each  packet  passing 
through  this  control  node  is  marked;  if  a  packet  tries  to  pass  through 
this  control  node  a  second  time,  it  is  typically  destroyed,  creating  an 
enpty  slot. 


3.4.3  DLCN.  The  link  level  transmission  scheme  for  the  distri¬ 
buted  loop  computer  network  (PLCN)  [6,8,9]  uses  a  shift  register  insertion 
technique  to  place  variable  (but  hardware  restricted)  length  messages 
onto  a  ring.  Two  shift  register  buffers  are  used;  one  is  a  variable 
length  delay  buffer  that  receives  data  from  the  predecessor  node,  and 
the  other  is  a  fixed  length  shift  register  that  contains  data  to  be 
placed  onto  the  ring  at  the  present  node  interface. 

A  message  arriving  for  transmission  at  a  given  node  waits  in  the 
output  buffer  until  end  of  message  is  detected  for  the  data  message  pass¬ 
ing  through  that  node  from  predecessor  to  successor  nodes.  When  this 
event  occurs,  new  incoming  data  from  the  predecessor  node  is  routed  into 
the  delay  buffer,  and  data  in  the  output  buffer  is  shifted  out  onto  the 
ring,  thereby  splicing  the  waiting  message  at  this  node  between  two  mes¬ 
sages  already  in  transit  on  the  ring.  In  other  words,  so  long  as  there 
is  enough  space  available  in  the  delay  buffer  to  hold  an  incoming  message, 
precedence  is  usually  given  to  transmitting  a  newly  arrived  or  already 
waiting  message  at  the  present  node  ahead  of  an  incoming  message  already 
on  the  ring.  This  technique  tends  to  minimize  waiting  times  for  messages 
to  be  placed  onto  the  ring  at  the  expense  of  randomly  delaying  transmitted 
messages  en  route  to  their  destinations.  The  maximum  length  message,  which 
is  in  effect  a  variable  but  maximum  length  packet,  is  fixed  by  the  length 
of  the  delay  buffer  at  each  node.  When  a  message  reaches  its  destination 
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It  Is  removed  from  the  ring  by  that  node.  If  the  message  Is  received 
correctly,  a  high  priority  acknowledgement  message  1$  placed  on  the  ring 
by  the  destination  node,  addressed  to  the  source  of  the  received  message. 

Presumably,  a  message  whose  source  or  destination  fields  are 
corrupted  will  be  error  checked  In  such  a  fashion  as  to  prevent  the 
wrong  destination  from  acknowledging  correct  receipt  of  the  message.  As 
with  receipt  of  a  negative  acknowledgement,  lack  of  a  positive  acknowledge¬ 
ment  after  some  appropriate  time  period  (called  a  time  out)  could  cause 
the  source  node  to  retransmit  the  data  message.  A  message  unclaimed  by 
Its  destination  would  also  presumably  be  removed  from  the  ring  when  the 
source  address  Is  recognized  by  some  source  node  as  part  of  Its  check 
and  forward  operations.  Since  DLCN  uses  a  distributed  control  mechanism, 
there  Is  no  central  controller  to  perform  any  of  these  functions. 

3.4.4  Playthrough.  The  Playthrough  mechanism  for  distributed 
control  of  r)ng  networks  [10,11,121  Is  a  check  and  forward  link  level 
control  protocol  that  provides  for  simultaneous  transmission  of  multiple 
variable  length  messages  of  any  length.  Control  Is  completely  distributed, 
and  data  and  control  messages  both  share  the  ring.  Control  Is  based  on 
a  special  synchronizing  message  (or  token)  called  GO  that  differs  from 
the  Newhall  synchronizing  token  In  two  ways:  first,  60  precedes  rather 
than  follows  data  messages,  so  that  It  can  continue  around  the  ring 
seeking  new  messages  to  activate;  and  second,  GO  circulates  perpetually 
despite  the  presence  of  other  traffic.  This  perpetual  circulation  Is 
achieved  by  giving  GO  a  higher  priority  and  allowing  It  to  preempt  tenpo- 
rarlly  any  data  message  It  overtakes.  Thus  GO  appears  at  times  to 
travel  Inside  data  messages,  or  In  golfing  terms,  to  "playthrough."  The 
protocol  bears  the  name  of  this  distinctive  feature. 

When  GO  arrives  at  any  node  with  a  message  to  send,  transmission 
may  begin  If  there  Is  a  free  path  to  the  destination.  To  Implement  this 
rule  without  collisions,  other  control  messages  precede  and  follow  the 
data  message  to  update  the  other  nodes  about  changes  In  loop  status.  Thus 
the  nodes  trust  be  able  to  recognize  control  messages  and  maintain  a  modest 
amount  of  local  Information  about  the  ring.  In  order  to  propagate  such 
status  Information,  the  update  control  messages  play  through  any  data 
messages  they  encounter.  Although  the  update  messages  are  synchronized 
by  GO,  their  even  higher  priority  causes  them  to  precede  GO  so  that  each 
node  has  the  correct  status  Information  before  GO  arrives. 

Some  operational  aspects  of  this  ring  are  worth  noting.  Data 
messages  can  be  preempted  only  at  their  sources.  This  means  that  there 
Is  no  store  and  forward  phenomenon  or  buffering  delay  en  route  to  the 
destination,  except  for  a  small  fixed  amount  at  each  node.  The  delays 
from  preenptlon  are  brief  because  the  Intervening  control  messages  are 
short.  Hence,  the  primary  message  delay  Is  due  to  queueing  at  the  source. 

Except  for  GO  >/h1ch  continues  traveling,  each  control  message 
makes  exactly  one  complete  circuit  of  the  ring  and  Is  removed  by  Its 
source.  This  permits  acknowledgements  from  the  destination  node  to  ride 
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for  free  on  returning  control  messages  and  to  avoid  queueing  delays.  In 
addition,  control  messages  complete  the  round  trip  in  a  fixed  time  that 
can  be  determined  dynamically.  This  enables  a  very  accurate  timeout  mecha¬ 
nism  to  be  used  for  error  detection  and  for  capture  and  removal  of  unacknowl 
edged  or  corrupted  control  messages. 

3.5  GPSS  Models  of  Ring  Networks.  Program  Modifications,  and  Correc- 
tions. 


Three  network  models  written  by  C.C.  Reames  [6]  in  GPSS/360  were 
obtained  through  the  assistance  of  Professor  M.T.  Liu  of  The  Ohio  State 
University.  These  programs  for  the  Newhall,  Pierce,  and  DLCN  single  ring 
conputer  networks  were  then  modified  to  run  under  GPSS/360  on  the  APG  IBM 
360/65.  Listings  of  these  GPSS/360  programs  can  be  found  in  the  appendices 
of  the  PhD  dissertation  by  ReamesC6],  pages  178  to  194.  Short  excerpts 
showing  our  modifications  to  these  programs  are  included  in  Appendix  B. 

They  were  also  translated  into  GPSS  1100  for  execution  on  the  ARRAOCOM 
UNIVAC  1108.  GPSS  1100  listings  can  be  found  in  Appendix  C;  the  line  for 
line  comments  are  the  same  as  those  for  the  IBM  versions  in  [6]  and  were 
thus  omitted  here. 

A  GPSS  1100  simulation  program  appearing  in  [7]  for  the  Play- 
through  protocol  ring  network,  found  here  in  Appendix  C,  was  modified  and 
corrected  slightly  and  also  translated  into  the  GPSS/360  version  found  in 
Appendix  B.  In  this  case,  line  for  line  comments  are  included  in  the  GPSS 
1100  and  the  GPSS/360  versions  to  align  the  translations. 

Several  modifications  to  the  original  6PSS/360  and  GPSS  1100 
programs  were  made;  some  changes  were  necessary  to  allow  the  models  to 
execute  under  GPSS/360  and/or  GPSS  1100,  and  some  were  made  to  align  the 
assumptions  concerning  message  routing  and  error  handling  and  to  correct 
minor  errors. 


“^.S.l  Changes  to  the  Pierce  Model.  The  GPSS/360  (enhanced) 

Pierce  network  simulation  program  referred  to  the  absolute  clock  standard 
numerical  attribute,  which  is  not  available  directly  in  either  of  the 
available  versions  of  GPSS/360  or  GPSS  1100.  Hence,  additional  code  to 
effectively  simulate  the  absolute  clock  facility  was  placed  into  the  Pierce 
network  simulation  programs  between  labels  LASTP  and  PATW. 

3.5.2  Changes  to  the  DLCN  Model.  The  original  DLCN  program  [6] 
attenpts  to  simulate  the  effects  on  system  loading  and  total  message  transit 
time  (or  end  to  end  delay)  caused  by  noise  corrupted  messages  that  include 
one  or  more  erroneous  characters.  If  the  message  (i.e.  transaction)  is 
marked  as  being  received  in  error,  it  is  discarded  by  the  destination 
node,  a  negative  acknowledgement  is  sent  to  the  source  node,  and  the 
message  is  placed  at  the  front  of  the  source  node  message  queue  for  retrans¬ 
mission.  Unfortunately,  the  implementation  of  this  feature  incorrectly 
counts  the  erroneous  message  as  a  successful  reception  (in  terms  of  the 
statistics  for  end  to  end  delay,  and  queueing  time),  and  then  resets  the 
corresponding  message's  time  in  system  to  zero  so  that  it  appears  and  is 
counted  in  the  statistics  as  a  newly  arriving  message  that  encounters 
hardly  any  queueing  time  thereby  slightly  skewing  the  output  statistics. 
Because  of  this  approach  to  handling  the  simulation  Of  erroneous  messages 
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with  a  mean  character  error  rate  of  one  In  ten  thousand,  mean  total  trans¬ 
mission  time  for  all  messages  handled  ty  the  network  when  errors  are  per¬ 
mitted  to  occur  Is  about  10  percent  lower  than  the  mean  total  transmission 
time  found  when  no  errors  occur,  as  seen  In  Figure  3.2.  Such  a  result  Is 
counterintuitive  and  slightly  Incorrect.  Because  the  other  ring  network 
simulation  progr'^ms  have  no  provisions  for  handling  messages  with  errors 
In  them,  the  character  error  generation  facility  in  the  DLCN  program  was 
disabled,  resulting  In  a  version  referred  to  as  DLCNNE  for  "no  errors". 

This  allows  a  more  uniform  comparison  of  simulation  results  for  the  dif¬ 
ferent  ring  network  protocols  and  removes  an  apparent  cause  of  skewed 
results  in  the  total  time  statistics  for  DLCN. 

3.5.3  Changes  to  the  Playthrough  Model.  Because  of  the  rather 
complex  and  specific  ordering  in  vd^lch  messages  must  be  placed  on  the 
communication  links,  the  Playthrough  simulation  program  maintains  Its  own 
user  chains,  which  are  in  effect  user  controlled  transaction  queues.  The 
user  chain  is  scanned  In  first-in  first-out  order  to  locate  the  first 
message  in  the  queue  having  a  free  path  to  its  destination.  If  one  is 
found,  that  message  (transaction)  Is  removed  from  the  queue  and  the  remain¬ 
ing  entries  are  left  In  their  original  order  in  the  queue.  This  is  accom¬ 
plished  by  circularly  shifting  the  queue  entries  and  examining  the  leading 
entry  until  either  a  message  with  a  free  path  is  found  or  until  the  queue 
is  restored  to  its  original  condition  given  the  number  of  elements  on  the 
queue.  Sobel's  original  Implementation  for  certain  queue  conditions  mis¬ 
counted  the  number  of  circular  shifts  by  one  so  that  reordering  of  the 
queue  after  removal  of  an  Interior  entry  left  one  element  out  of  position, 
resulting  in  occasionally  increased  waiting  times  for  some  transactions. 

A  minor  modification  to  the  logic  governing  chain  reordering  corrected 
this  problem. 

The  Playthrough  message  destination  assignment  scheme  was  modi¬ 
fied  to  match  that  found  in  the  Reames  models  so  that  the  distribution  of 
destinations  is  uniform.  Sobel's  original  scheme  generated  message 
destinations  skewed  toward  shorter  distances. 

3.6  GPSS  Ring  Network  Simulation  Results. 

This  section  describes  the  results  of  running  both  IBM  GPSS/360 
and  UNIVAC  GPSS  1100  programs  for  the  various  ring  network  models.  It  was 
assumed  that  published  data  [9]  were  based  on  the  same  startup  and  run 
termination  conditions  found  In  the  Reames  programs  from  [6].  The  Pierce 
model  uses  a  startup  of  250  messages  to  preload  the  queues  and  initialize 
the  system,  and  then  accumulates  statistics  on  the  successful  transmission 
of  1200  additional  messages.  The  Newhall  model  uses  a  startup  of  200  mes¬ 
sages  for  Initialization  and  then  accumulates  statistics  over  1000  addi¬ 
tional  messages.  The  DLCN  and  Playthrough  models  both  use  a  startup  of 
100  messages  and  accumulate  statistics  on  1000  successfully  transmitted 
messages.  Without  detailed  statistical  analyses  of  these  startup  para¬ 
meters  to  determine  if  steady  state  has  actually  been  reached,  these 
seemingly  arbitrary  but  Intuitively  justifiable  choices  lead  to  acceptable 
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FOR  SIX-NODE  DLCN  LOOP  NETWORK  WITH 
NOMINAL  MEAN  MESSAGE  LENGTH  OF 
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Figure  3.2.  Comparison  of  DLCN  Output  Data  with  those 
for  DLCNNE. 


qualitative  results  only  If  one  Is  Interested  In  gaining  an  Idea  of  rela¬ 
tive  performance  differences.  A  check  of  startup  conditions  plotting 
relative  changes  In  the  mean  of  output  parameters  was  made  for  DLCN  and 
Playthrough  Indicating  that  100  terminated  messages  seems  to  be  sufficient 
for  the  warmup  period.  However,  If  one  wishes  to  draw  statistically 
valid  Inferences  from  the  simulation  results,  one  should  use  formal 
statistical  tests  while  collecting  the  data. 

Three  Inportant  areas  of  concern  are  (1)  starting  criteria  for 
data  collection,  (2)  stopping  criteria,  and  (3)  determining  to  what  degree 
the  data  are  correlated.  Starting  criteria  are  concerned  mainly  with 
determining  at  what  point  the  simulation  closely  approximates  steady- 
state.  Stopping  criteria  determine  when  (how  soon)  It  Is  statistically 
safe  to  stop  collecting  data  and  still  be  able  to  draw  conclusions  with 
the  required  level  of  confidence.  Correlated  data  yield  less  Information 
about  a  system  per  observation  than  If  all  data  were  Independent.  To 
compensate  for  this  lower  average  Informational  content,  one  must  collect 
more  data.  Later  simulations  of  DLCN  by  Itself  for  example  [211  take 
cognizance  of  these  Items.  Although  statistical  validity  of  simulation 
results  was  not  the  main  concern  of  this  study.  It  must  be  a  major 
consideration  of  any  production  oriented  simulation  study  on  whose  results 
decisions  are  to  be  based. 

3.6.1  Message  Interarrival  Time  and  Length  Distributions. 

Tests  of  tbe  correctness  of  generated  exponential  olstrlbutlons  In  both 
IBM  and  UNIVAC  simulations  were  performed.  Because  message  arrivals  at 
each  netv/ork  node  (from  Its  attached  component)  are  assumed  to  be  governed 
by  a  Poisson  process  with  Identical  parameters  at  each  node,  plots  of 
actual  Interarrival  times  were  made  to  see  If  they  resemble  exponential 
distributions  and  to  see  If  those  generated  by  the  UNIVAC  Intrinsic 
exponential  function  are  similar  to  the  IBM  user  defined  exponential 
function.  One  such  example  plot  showing  count  of  the  number  of  messages 
versus  corresponding  Interarrival  time,  where  Interarrival  times  are 
grouped  Into  ten  unit  Intervals,  Is  shown  In  Figure  3.3.  The  mean  Inter¬ 
arrival  time  Is  300  character  times  at  each  node;  for  the  six  node  system 
considered  here  the  system's  mean  Interarrival  time  Is  300  divided  by  6. 

A  sample  plot  of  the  count  of  the  number  of  messages  versus 
corresponding  message  length  Is  shown  In  Figure  3.4.  Each  generated 
message  has  nine  characters  of  overhead  Information  added  to  Its  length, 
and  each  frequency  count  was  accumulated  over  a  ten  character  Interval 
after  overhead  Information  was  appended;  hence,  the  first  Interval  counts 
messages  of  length  between  nine  and  ten  characters  only  thus  skewing  the 
plot  from  a  true  exponential.  All  of  the  ring  network  simulation  programs 
considered  here  use  an  approximate  exponential  distribution  for  generating 
message  lengths  truncated  at  a  maximum  of  500  characters  because  of  the 
DLCN  hardware  defined  delay  buffer  limit  of  512  characters  Including 
overhead. 

Overall,  the  UNIVAC  and  IBM  generators  produce  similar  results 
for  exponentially  distributed  Interarrival  times  and  message  lengths. 
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Figure  3.3.  Generated  Arrival  Distribution  (DLCN). 


MLEN  =  50 


31 


Figure  3.4.  Message  Length  Distnbution  (DLCN). 


The  IBM  plots  are  based  on  a  sixty  point  user  defined  continuous  approxi¬ 
mation  function,  and  the  IJNIVAC  plots  are  based  on  the  GPSS  1100  intrinsic 
exponential  function. 

3.6.2  Some  Effects  of  Varying  Pseudo  Random  Number  Sequences. 
Plots  of  end  to  end  delay  (or  total  transmission  time)  versus  messaoe 
arrival  rate  at  each  node  are  shown  in  Fiqure  3.5;  the  two  plots  shoivn  are 
for  both  IP.M  and  UMIVAO  sirwlations  of  the  six  node  OLCN  ring  using 
different  combinations  of  pseudo  random  number  generators  (or  the  same 
generator  in  the  IBN!  case  differently  seeded).  For  runs  of  1000  message 
terminations,  the  end  to  end  delay  is  obviously  sensitive  to  the  seguences 
of  pseudo  random  numbers  used.  To  smooth  these  differences  one  can  make 
several  simulation  runs  using  either  different  sets  of  random  number 
generators  or  different  sets  of  seeds  and  then  either  take  the  mean  of  the 
results  associated  with  each  designated  interarrival  time,  or  construct 
the  final  curve  using  minimum  mean  square  error  fit.  This  would  be  the 
case  if  fixed  termination  counts  are  used  or  if  the  simulation  is  stopped 
at  a  fixed  time.  A  statistically  better  approach  v/ould  be  to  design  the 
stopping  criteria  to  take  cognizance  of  the  confidence  intervals  involved 
with  the  statistics  of  the  output  data,  as  mentioned  earlier. 

3.6.3  Nominal  Versus  Measured  Parameters.  Differences  were 
observed  in  UNIVAC  and  IPM  flPSS  outputs  for  DLFN  simulations  using  identi¬ 
cal  nominal  parameters  for  both  mean  message  length  and  mean  interarrival 
time.  The  differences  in  observed  mean  message  lengths  are  essentially 
constant  for  all  corresponding  interarrival  times  (for  UfllVAC  a  mean  of 
5B.4^  0.2  characters  and  for  IBM  a  tnean  of  57.9  ^  0.1,  making  the  v/orst 
case  difference  approximately  It  of  nominal  mean  of  59  characters  includ¬ 
ing  the  9  character  overhead).  Because  the  differences  in  observed  mean 
message  length  are  essentially  constant,  only  differences  in  mean  interar¬ 
rival  times  appear  to  be  significant.  For  the  six  node  OLCN  simulation 
with  a  nominal  mean  message  length  of  50  characters  (excluding  overhead) 
t\Ki  curves  are  shov/n  in  Figure  3.6  for  both  IBM  and  UNIVAC  simulation 
results  for  total  message  transmission  time  (i.e.,  end  to  end  delay). 

The  curves  marked  "nominal"  are  plotted  using  the  nominally  specified 
nodal  interarrival  times.  The  curves  marked  "adjusted"  use  an  abscissa 
of  observed  mean  nodal  interarrival  times.  The  "total"  time  ordinates 
using  nominal  interarrival  time  values  are  skev/ed  to  the  high  side  for 
the  UNIVAC  results  and  are  skev/ed  slightly  to  the  lo»y  side  for  the  IBM 
results,  thereby  giving  a  more  pessimistic  estimate  of  system  performance 
for  UNIVAC  data  and  a  more  optimistic  estimate  of  perforiTiance  for  IBM 
data  than  is  the  case  if  observed  mean  interarrival  times  are  used  as 
abscissas. 


3.6.4  Results  for  Newhall  Loop.  Simulation  results  for  total 
transmission  times  versus  per  node  message  arrival  rate  for  the  Nev/hall 
Loop  Network  are  shown  in  Figure  3.7.  Curves  for  IBM  GPSS/36n,  UNIVAC 
GPSS  1100,  and  the  published  data  of  Peames  and  Liu  [9’’  are  shown  for 
comparison.  The  differences  are  likely  caused  by  variations  in  the 
actual  pseudo  random  number  sequences  used  in  each  case  coupled  with  the 
1000  transmitted  messages  stopping  criterion.  The  ir>M/360  and  UNIVAC 


32 


results  match  reasonably  well,  indicating  that  successful  and  correct 
translation  between  syntactically  different  dialects  of  GPSS  is  feasible. 

3.6.5  Comparison  of  Results  For  All  Four  Networks.  Paralleling 
the  stuCly  reported  In  L9!,  the  primary  quantities  of  interest  in  this 
study  are  the  mean  total  transmission  time  for  messages  (i.e.,  end  to 
end  delay)  and  mean  queueing  time  for  messages;  however,  many  other 
quantities  such  as  communication  link  utilizations  were  also  measured 
in  these  sinulations.  Some  of  the  relevant  times  that  are  discussed 
further  are  defined  below: 

(1)  queueing  time--time  elapsed  from  message  generation  until 
placement  on  the  loop  by  the  transmitter  at  the  source  node; 

(2)  transmission  t1me--time  elapsed  from  message  placement  on 
the  loop  until  the  last  character  is  received  and  removed  from  the  loop 
at  the  destination; 

(3)  acknowledgement  time — time  elapsed  from  generation  of  the 
acknowledgement  message  at  the  destination  node  until  the  last  character 
is  received  at  the  source  node; 

(4)  total  message  transmission  time  (or  end  to  end  delay  time)-- 
sum  of  (1)  and  (2)  only  for  Newhall  and  Pierce  loops;  sum  of  (1),  (2) 

and  (3)  for  DLCN  (including  DLCNNE);  and  modified  sum  of  (1),  (2)  and 
(3)  for  Playthrough,  where  Playthrough's  simulation  differs  from  the 
others  in  that  detailed  simulation  of  character  by  character  transmission 
does  not  take  place;  rather,  [control  message— data  message— control 
message!  groupings  of  characters  are  used  for  simulator  efficiency,  and 
the  acknowledgement  rides  for  free  on  the  trailing  control  message. 

(Note  that  inclusion  of  character  error  simulations,  not  currently  used 
in  any  of  the  loop  network  simulations,  would  likely  require  modifi¬ 
cation  of  Playthrough  code  to  perform  character  by  character  trans¬ 
mission  between  transmitter-receiver  pairs  around  the  ring  in  a  fashion 
similar  to  the  other  three  simulation  models.  Such  modification  would 
also  tend  to  increase  the  running  time  for  the  Playthrough  simulation.) 

The  general  characteristics  of  all  four  networks  modeled  are 
the  same.  Each  comprises  six  nodes,  with  each  message  source  being  an 
identical  and  independently  distributed  Poisson  process.  Messages  produced 
at  each  node  are  addressed  uniformly  to  the  other  five  nodes,  so  that 
message  traffic  is  entirely  symmetric  and  random.  Message  data  lengths 
are  assumed  to  be  exponentially  distributed  with  a  nominal  mean  of  50 
characters;  actually,  a  truncated  exponential  distribution  is  used  with 
no  message  exceeding  500  characters  in  length  in  order  not  to  violate 
the  hardware  defined  maximum  length  message  including  overhead  of  512 
characters  that  DLCN  was  assumed  capable  of  handling.  For  the  three 
loops  other  than  Playthrough,  nine  additional  characters  of  header  infor¬ 
mation  were  added  to  each  message  or  packet  produced;  the  Playthrough 
simulation  adds  ten  characters  of  overhead  in  the  following  way:  three 
characters  of  control  message  information  to  initiate  transmission  on 
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the  loopt  four  characters  of  overhead  added  to  the  data  message  In  the 
form  of  two  characters  of  message  length  Information  and  two  characters 
for  error  detection,  finally  followed  by  three  characters  of  control 
Information  to  terminate  the  loop  connection  from  source  to  destination 
and  to  carry  acknowledgement  Information  from  destination  to  source. 

All  timing  Is  In  arbitrary  character-time  units,  so  that  no  particular 
line  rate  Is  assumed.  Propagation  delay  on  the  conminicatlon  channel 
itself  was  Ignored.  In  the  three  models  other  than  Playthrough  each 
ring  Interface  unit  through  which  messages  pass  contributes  two  units  of 
delay:  one  unit  In  the  receiver  for  address  checking  and  one  unit  In 
the  transmitter.  In  Playthrough  GO  Is  delayed  by  only  one  time  unit  In 
ring  interfaces  with  nothing  to  transmit,  and  Is  delayed  by  three  time 
units  when  preceded  by  a  three  character  control  message  to  allow  time 
for  address  checking  and  control  message  transformation  at  appropriately 
designated  nodes  before  relay  by  the  ring  interface  transmitter.  Special 
features  in  the  DLCN  model  are  described  further  in  [9]. 

Tables  3.1  through  3.4  present  relevant  results  of  the  simula¬ 
tions  for  the  four  ring  networks  under  consideration.  In  all  of  these 
tables  certain  abbreviations  are  common  and  are  discussed  in  this  paragraph 
'tore  specific  labels  and  names  relating  to  measured  quantities  and  names 
used  in  the  program  listings  in  Appendices  B  and  C  are  discussed  in 
corresponding  specific  subparagraphs  below.  The  first  column  In  each 
table  lists  the  nominal  mean  message  interarrival  time  at  each  node  in 
the  corresponding  network.  Again,  the  units  are  character-times. 

Because  the  network  (or  system)  comprises  six  nodes,  the  nominal  mean 
system  interarrival  time  is  one  sixth  of  this  value.  The  third  and 
fourth  columns  display  both  mean  and  standard  deviation  of  the  measured 
system  interarrival  times  as  tabulated  in  the  programs  using  the  symbolic 
name  MSGAR  In  the  Newhall,  Pierce,  and  Playthrough  models,  and  the  name 
GENAR  in  the  OLCNNE  model.  The  node  arrival  rate  shown  In  column  two  Is 
confuted  as  the  reciprocal  of  six  times  the  mean  system  interarrival  time 
value  from  column  three.  Columns  five  and  six  show  means  and  standard 
deviations  for  the  measured  mean  message  lengths  (with  program  name 
MSGLN).  The  reasons  these  values  differ  significantly  from  the  nominal 
mean  of  fifty  characters  are  due  to  both  underestimation  of  target  mean 
by  the  truncated  exponential  using  the  IBM  pseudo  random  number  generator 
and  to  the  way  in  which  header  characters  are  accounted  for,  as  discussed 
in  the  model  specific  paragraphs  below.  The  seventh  column  lists  mean  fad 
lity  utilization  which  is  found  by  averaging  the  six  facility  mean  utiliz¬ 
ations.  Each  facility  (or  transmitter)  utilization  essentially  measures 
utilization  of  the  corresponding  outgoing  communication  link. 

3. 6. 5.1  Model  Specific  Items  for  Newhall.  Table  3.1  displays 
means  and  standard  deviations  of  two  simulation  output  parameters  of 
Intense  Interest  and  a  third  of  only  moderate  Interest.  The  mean  total 
queueing  time  experienced  by  all  messages  arriving  for  transmission 
anywhere  In  the  system  of  six  nodes  is  tabulated  in  the  simulation  model 
under  the  name  TLQTM  and  is  listed  in  Table  3.1  as  one  entry  for  each 
corresponding  message  interarrival  time.  Message  transmit  time  Is  shown 
under  the  heading  TRNTM,  and  total  message  transmission  time  which  is 


approximately  the  sum  of  TLQTM  and  TRNTM  (though  tabulated  separately 
in  the  model)  is  shown  under  the  heading  TMGTM.  The  measured  mean  message 
length  tabulated  under  heading  MSGLN  is  based  on  a  nominal  mean  message 
length  of  50  plus  9  overhead  characters  (or  59  characters). 

3. 6. 5. 2  Model  Specific  Items  for  Pierce.  For  the  Pierce  loop 
simulation  results  the  measured  mean  message  length  is  nearly  the  nominal 
mean  value  of  50  characters.  The  nine  character  overhead  is  added  to 
each  packet  which  consists  of  at  most  36  characters,  and  the  average 
number  packets  per  message  (NPKMG)  is  2.35.  The  average  packet  synchro¬ 
nization  time  (SYNTM)  is  17.4;  the  average  packet  transmit  time  (PTRTM) 

is  46.6,  with  standard  deviations  shown  in  Table  3.2.  The  columns  labeled 
PKWTM  display  packet  waiting  time  statistics,  and  under  TPKTM  display 
total  packet  transmit  time.  The  main  parameter  of  interest  is  the  total 
message  transmission  time  displayed  under  TMGTM. 

3.6.5.3  Model  Specific  Items  for  DLCNNE.  Measured  mean  message 
length  for  the  PLCN  simulation  with  character  error  generation  facilities 
disabled  as  shown  in  Table  3.3  is  based  on  a  nominal  mean  length  of  50 
characters  plus  nine  characters  of  overhead.  Means  and  standard  deviations 
for  the  following  parameters  of  interest  are  displayed  in  the  remaining 
columns  of  Table  3.3.  Statistics  for  total  queueing  time  are  shown 

under  heading  TRQTM;  those  for  total  transmit  time  for  data  messages  on 
the  way  to  their  destinations  is  shovrt?  under  RCVTM,  and  total  transmit 
time  for  the  return  acknowledgement  message  is  shown  under  ACKTM.  TLATM 
is  the  total  message  transmission  time  which  is  (approximately)  the  sum 
of  TRQTM,  RCVTM,  and  ACKTM,  and  it  is  this  value  that  is  plotted  in 
Figure  3.8.  DLYTM  records  statistics  for  the  per  node  time  messages 
spend  in  delay  buffers  enroute  to  their  destinations. 

3. 6.5. 4  Model  Specific  Items  for  PI ayth rough.  Measured  mean 
message  lengths  shov/n  in  Table  3.4  for  the  Playthrougn  model  are  based 
on  a  nominal  mean  message  length  of  50  plus  4  overhead  characters  for  a 
total  of  54  characters.  The  six  additional  control  message  characters 
needed  to  start  and  stop  data  message  transmissions  affect  queueing  and 
total  time  statistics,  but  are  not  included  in  the  message  length  statis¬ 
tics.  Only  the  parameters  of  greatest  interest  are  shown  in  Table  3.4, 
namely  total  queueing  time  under  heading  TLQTM  and  total  transmission 
time  (plotted  in  Figure  3.8)  shown  under  heading  TTLTM.  TTLTM  includes 
the  acknowledgement  time  embedded  in  the  control  mechanism.  (Note, 
message  transit  time  is  the  difference:  TTLTM  minus  TLQTM.)  Table  3.5 
displays  additional  information  for  the  Playthrough  loop,  whe’^e  mean 
queueing  times  versus  distance  (in  number  of  nodes  to  the  destination) 
are  tabulated.  Average  waiting  times  for  messages  with  destinations  one 
hop  away  are  shown  under  heading  TLQl,  and  those  for  messages  with  desti¬ 
nations  five  nodes  away  are  shown  under  heading  TLQ5.  The  maximum  number 
of  messages  waiting  in  any  of  the  six  queues  as  well  as  the  average  number 
of  messages  waiting  in  queue  during  the  simulation  are  also  tabulated 
against  corresponding  message  interarrival  times  per  node. 
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PIERCE  LOOP  SIMULATION  RESULTS 
NOMINAL  MEAN  MESSAGE  LENGTH  =  50  CHARACTERS 
PACKET  LENGTH  =  36  CHARACTERS 


OLCNNE  (NO  TRANSMISSION  ERRORS)  LOOP  SIMULATION  RESULTS 
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PLAYTHROUGH  LOOP  SIMULATION  RESULTS 
NOMIMAL  MEAN  MESSAGE  LENGTH  =  50  CHARACTERS 
TIME  UNITS  ARE  "CHARACTER  TIMES" 


TABLE  3.5 


203.8  657.1  1982.7  3620.7  4924.8  27  8.340 


B 


3.7  Findings. 

The  data  generated  for  Newhall,  Pierce  and  DLCNNE  loops  agree 
reasonably  well  with  published  data  [91  in  that  the  relative  positions  of 
the  plotted  total  transmission  time  data  are  similar.  The  exact  values 
differ  somewhat,  which  for  Newhall  and  Pierce  can  be  accounted  for  by 
pseudo  random  number  generator  variations.  DLCNNE  differs  from  DLCN 
results  because  of  the  disabling  of  the  erroneous  message  generation  and 
retransmission  scheme  resulting  in  an  approximately  ten  per  cent  difference 
in  computed  values  as  discussed  in  Section  3.5.2  of  this  report. 

The  significance  of  Figure  3.8  is  that  it  provides  the  first 
extensive  comparison  between  the  DLCN  and  Playthrough  link  level  protocol 
schemes.  Overall  transmission  times  for  DLCN  are  lov/er  on  the  average 
than  for  the  other  link  level  protocol  schemes,  and  this  is  to  be  ex¬ 
pected.  tinder  heavy  loading,  the  Newhall,  Playthrough,  and  even  Pierce 
schemes  suffer  from  increased  queueing  delays,  whereas  the  DLCN  scheme 
is  designed  to  minimi?e  queueing  delays.  Nothing  is  free,  hov/ever,and 
in  the  DLCN  scheme  messages  suffer  random  exponentially  increasing  de¬ 
lays  en  route  to  their  destinations,  making  strict  timeouts  for  error 
control  difficult.  Transit  times  in  Playthrough  grow  approximately 
linearly  with  almost  imperceptible  slope,  so  that  as  in  Newhall,  once  a 
message  transmission  is  initiated  it  proceeds  rapidly  and  is  completed 
in  almost  fixed  time.  The  disadvantage  of  Playthrough  is  that  under 
heavy  load,  queueing  time  gro\-/s  exponentially  because  long  hop  messages 
must  wait  long  times  before  a  sufficient  number  of  links  from  source  to 
destination  nodes  become  simultaneously  free. 

These  disadvantages  are  common  among  schemes  that  use  dedicated 
circuit  switching  in  the  transmission  of  messages.  Packet  switched 
schemes  tend  to  experience  less  rapid  growth  in  queueing  time  under  heavy 
loads;  hov/ever,  they  require  dedicated  intelligence  or  capacity  in  either 
the  ring  interface  processor  or  in  the  attached  component  (e.g. ,  the  host 
computer)  to  packeti/e  messages  at  their  sources  and  to  reassemble  at  their 
destinations  packets  that  are  arriving  in  arbitrary  sequence  from  possibly 
disparate  messages.  DLCN  employs  variable  length  packets  in  this  sim¬ 
ulation  up  to  a  maximum  of  512  characters  in  length,  which  represents  a 
chosen  hardiyare  limit.  Messages  of  longer  length  were  not  allowed  in  this 
simulation  because  the  code  to  packetize  them  was  not  included  in  the 
model.  DLCN  minimizes  queueing  times  by  usually  placing  on  the  ring 
nev/ly  arriving  messages  ahead  of  messages  already  on  the  ring  through  the 
use  of  expandable  delay  buffers.  This  technique  appears  to  be  a  particu¬ 
larly  effective  means  of  maintaining  reasonable  mean  transmission  times 
under  heavier  loads  than  is  possible  with  the  other  link  level  schemes. 

An  advantage  of  the  Nev/hall  and  Playthrough  protocols  is  their  ability  to 
transmit  quickly  messages  of  any  arbitrary  length  ivhen  the  number  of 
characters  arriving  for  transmission  to  the  entire  network  does  not 
exceed  the  burst  character  transmission  rate. 

An  interesting  observation  from  examining  the  plots  in  Figure 
3.8  is  that  the  perpetually  circulating  control  token  in  the  Playthrough 
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schene  tends  to  have  a  packetizino  effect  on  mean  total  transmission 
tines;  so  that  for  non-saturatinq  loads  it  corresponds  to  hut  is  lower 
than  the  mean  total  transmission  tine  for  the  Pierce  schene. 

Neither  the  Pierce  nor  the  Hewhall  simulations  include  the  load- 
inq  effects  and  delays  produced  hy  the  inclusion  of  acknowledqenents  for 
messaqes  sent;  whereas,  Playthrouoh  and  niCNNE  do  include  then.  The 
simulation  results  are  therefore  favorably  biased  for  Pierce  and  Mewball. 

4.  CONCLUSIONS  AHO  PECOHMENOATIOHS  FOR  FURTHER  WORK 

4.1  GPSS  Capabilities. 

Some  capabilities  of  the  GPSS  lanquaqe  for  model inq  and  sinu- 
latinq  systems  were  presented  in  Chapter  2,  and  differences  in  two  avail¬ 
able  implementations  of  this  lanquaqe  were  discussed.  GPSS  has  several 
facilities  for  system  level  model inq  of  computer  communication  networks. 
Messaqes  are  easily  modeled  as  dynamic  entities,  called  transactions. 
Lanquaqe  features  are  provided  for  qeneratinq  or  inplementinq  messaqe 
arrivals  and  other  randomly  occurinq  events  such  as  link  and  node  failures 
and  dynamic  routinq  schene  choices.  Equipment  entities  such  as  trans¬ 
mitters,  receivers  and  messaqe  queues  are  easily  modeled.  Both  automatic 
and  user  specified  means  for  collectinq  and  comnutinq  statistics  for 
messaqe  transmission  such  as  mean,  variance,  and  distributions  (per¬ 
centiles)  of  queueinq,  transmission,  and  end  to  end  delay  times  are  also 
included.  These  statistics  can  be  used  to  predict  system  behavior  under 
varyinq  conditions. 

4.2  Sample  Simulation  Results  and  Applications. 

To  illustrate  use  of  these  capabilities,  GPSS  models  of 
several  rinq  tonolonv  computer  conmiini cation  networks  were  examined  in 
Chapter  3.  Oata  were  collected  to  indicate  performance  under  varyinq 
system  load  for  each  of  the  rinq  network  link  level  protocols  considered. 
These  oerformance  data  were  plotted  to  show  relative  performance  of  the 
differinq  link  control  and  messaoe  handlino  schemes.  Tests  for  statis¬ 
tical  validity  of  these  data  should  be  performed  before  decisive  con¬ 
clusions  are  drawn  from  these  comparisons.  It  was  the  purpose  of  this 
study  to  demonstrate  use  of  GPSS  for  computer  communication  network 
modelinq  rather  than  to  produce  statistically  valid  system  comparisons. 
However,  some  statistical  validity  tests  for  the  Playthrouqh  data  were 
performed  jisino  the  method  of  batch  means  employed  by  Wolf  [211  in  his 
simulation  of  a  double  loop  flLCN  configuration.  For  instance,  the  mean 
total  transmission  time  entries  (TTLTM)  in  Table  3.4  satisfy  a  RO* 
confidence  level  test  at  nominal  interarrival  times  of  300,  600,  and  1000. 
This  suqqests  reasonable  accuracy  in  the  plotted  performance  data. 

The  simulations  that  have  been  run  have  used  a  nominal  mean 
messaqe  lenqth  of  50  characters  (some  rwssaqes  are  longer  and  some  are 
shorter).  This  mean  messaqe  lenqth  approximates  the  characteristics 
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of  many  actual  computer  communication  schemes,  hut  the  actual  distributions 
involved  as  well  as  their  means  may  vary  some\^at  from  this  choice.  To 
gain  an  appreciation  of  how  well  the  ring  network  schemes  considered 
here  might  work  In  a  typical  computer  communications  structure,  we  must 
make  additional  assumptions  about  mean  interarrival  times  for  messages, 
the  number  of  binary  digits  or  bits  used  to  encode  the  characters, and 
the  speed  of  the  communication  links  in  the  network  In  bits  per  second 
to  make  it  independent  of  modulation  scheme. 

In  order  to  use  the  data  presented  in  Chapter  3  recall  we 
have  assumed  that  messages  are  on  the  average  50  and  no  more  than  500 
characters  in  length  with  length  governed  by  a  truncated  exponential 
distribution.  If  one  assumes  10  bits  are  required  to  transmit  one 
character  (7  code  hits,l  parity  bit,  1  start  bit,  and  1  stop  bit  for  an 
asynchronous  format),  then  one  “character  time"  at  a  link  transmitter/ 
receiver  speed  of  1  million  hits  per  second  (1  Mbps)  is  10-5  seconds, 
and  at  a  speed  of  1200  bps  is  8.33x10-3  seconds.  Assuming  a  netvmrk 
of  six  identical  data  terminals  ir  which  operators  send  messages  to  some 
destination  node  at  an  average  rate  of  one  every  30  seconds,  then  we 
compute  the  communications  netv/ork  (i.e.,  system)  arrival  rate  as:  multiply 
six  (the  number  of  nodes  corresponding  to  the  simulation  results  presented; 
times  the  per  node  arrival  rate  (in  messages  per  second)  times  the  time 
for  one  character  (in  seconds  per  character  time).  At  a  line  speed  of  1 
Mbps  these  assumptions  result  in  a  mean  system  arrival  rate  of  0.002x10-3 
messages  per  character  time  (or  0.00033  x  10-3  messages  per  character 
time  per  node).  Looking  in  Figure  3.8  under  this  arrival  rate,  one 
finds  that  for  all  four  ring  net\«rk  structures  considered  the  expected 
mean  total  transmission  time  for  messages  is  less  than  200  character 
times  which  corresponds  to  less  than  2  milliseconds.  At  a  link  speed 
of  12^0  bps  using  otherwise  same  assumptions,  the  per  node  arrival  rate 
is  2.78  X  10-3  messages  per  character  time.  At  this  arrival  rate 
Figure  3.8  says  that  for  all  but  the  Mewhall  scheme  the  expected  message 
transmission  time  is  less  than  540  character  times  (or  4.5  seconds);  for 
the  '‘lewhall  scheme  the  expected  message  transmission  time  is  approximately 
4000  character  times  nr  33  seconds,  not  a  very  desirable  performance 
if  one  expects  to  generate  a  new  message  for  transmission  once  every  30 
seconds. 


4.3  Simulation  Language  Alternatives. 

Having  the  capability  to  simulate  various  conputer  communication 
networks  quickly  permits  analysts  to  identify  potential  bottlenecks  and 
deficiencies  in  proposed  computer  communication  network  schemes.  Various 
discrete  event  languages  are  available  to  facilitate  the  programming  of 
these  simulation  models.  Two  of  the  more  popular  are  various  dialects 
of  PPSS  and  SPISCRIPT.  GPSS  is  a  block  oriented  language  in  which  simul¬ 
ator  specifications  relate  more  to  the  flow  of  dynamic  entities  in  the 
actual  model  than  to  traditional  computer  programming  languages.  GPSS 
is  interpreted  rather  than  compiled  as  is  SIMSCRIPT.  Various  comparisons 
of  those  languages  [23"!  and  [24’’  point  to  advantages  and  disadvantages 
of  each,  Reginners  usually  have  an  easier  time  learning  GPSS  because  of 
the  abundance  of  tutorial  material  available;  whereas,  far  less  complete 
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tutorial  material  is  available  to  beginners  learning  SIMSCRIPT.  Because 
SIMSCRIPT  is  compiled,  some  models  written  in  this  language  can  be  expec¬ 
ted  to  execute  more  rapidly  than  do  similar  models  written  in  GPSS,  Recent 
additions  to  the  SIMSCRIPT  language,  however,  tend  to  reduce  its  speed 
advantage  [23].  Certain  computations  are  more  easily  specified  in  one 
language  than  in  the  other.  For  instance,  exponentiation  is  not  available 
as  a  primitive  and  is  cumbersome  to  specify  in  GPSS  [21]  p.lll. 

4.4  Use  of  GPSS. 

Two  dialects  of  GPSS  (namely,  GPSS/360  and  GPSS  1100)  were  used 
in  the  example  computer  communication  network  simulations  documented  here. 
Differences  in  both  syntax  and  semantics  between  the  two  dialects  have 
been  identified  and  are  discussed  in  Chapter  2.  Because  of  these  differ¬ 
ences,  care  should  be  exercised  when  comparing  output  data  from  one  dialect 
with  that  from  another  in  order  to  insure  that  the  comparison  is  meaningful. 
It  is  however  possible  to  correctly  translate  a  model  from  one  dialect  to 
another  by  carefully  tracing  the  flow  of  transactions  in  the  two  models  to 
identify  and  correct  or  at  least  account  for  differences  in  interpreter 
execution  (i.e.,  semantics).  This  task,  however,  is  not  particularly 
easy  and  should  not  be  taken  lightly. 

Because  of  the  variety  of  programming  techniques  required  to 
implement  the  several  ring  network  simulation  models  in  GPSS,  the  collec¬ 
tion  of  programs  found  in  appendices  B  and  C  coupled  with  those  in  [6] 
should  be  a  valuable  aid  to  programmers  seeking  to  model  other  computer 
communication  network  architectures  and  protocols.  Each  protocol  considered 
has  its  own  peculiar  implementation  requirements  that  relate  to  other 
actual  and  potential  computer  network  structures. 

4.5  Future  Work. 

Because  of  recognized  deficiencies  in  the  GPSS  language,  such  as 
long  execution  times  and  cumbersome  constructions  to  do  simple  computa¬ 
tions  directly  available  in  other  languages,  an  investigation  into  the  use 
of  the  discrete  event  simulation  language  SIMSCRIPT  II. 5  should  be  considered 
for  further  work.  There  are  indications  that  SIMSCRIPT  II. 5  is  superior 
to  GPSS  because  of  its  generally  higher  speed  of  execution  and  lower  memory 
space  requirements  for  the  same  model  [23]  and  [24];  also,  model  implemen¬ 
tation  reportedly  requires  programmer  skill  roughly  equivalent  to  that  of 
a  competent  FORTRAN  or  ALGOL  programmer,  which  should  cause  little  diffi¬ 
culty  for  most  organizations.  The  set  of  examples  considered  in  Chapter  3 
could  be  used  as  a  starting  point(and  validation  check)  for  initial  SIMSCRIPT 

II. 5  modeling  efforts.  Use  of  SIMSCRIPT  will  not  necessarily  replace  the 
use  of  GPSS  because  some  investigators  [24]  indicate  that  it  is  likely  to 
be  faster  to  program  an  initial  system  model  in  GPSS  to  get  quick  results 
that  can  be  used  to  guide  the  development  of  a  more  comprehensive  (and 
possibly  more  efficient)  SIMSCRIPT  II. 5  model. 

Because  use  of  a  ring  network  architecture  has  been  proposed  for 
SIGMA  [29],  the  simulation  models  examined  here  should  be  considered  for 
potential  further  use  in  evaluation  of  the  SIGMA  computer  communications 
structure. 
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APPENDIX  A 


ON  THE  RANDOMNESS  OF  PSEUDO  RANDOM  NUMBER  GENERATORS 
USED  IN  IBM  GPSS/360  AND  UNIVAC  GPSS  1100  LANGUAGES 

A.l  INTRODUCTION 

Tests  of  randomness  were  performed  on  the  UNIVAC  GPSS  1100  and 
IBM  GPSS/360  pseudo  random  number  generators  when  simulation  models 
translated  from  one  language  to  the  other  failed  to  yield  conparable 
statistics  for  checkout  runs.  Initially,  the  translations  themselves 
were  suspect;  however,  subsequent  investigation  found  no  basis  for 
faulting  the  translations. 

The  simulation  models  tested  rely  on  pseudo  random  number  gener¬ 
ators  embedded  in  the  languages  to  generate  message  traffic  for  input  to 
the  models.  Small  differences  in  mean  message  lengths  and  mean  interar¬ 
rival  times  for  this  traffic  were  observed  for  corresponding  runs  in 
the  two  languages,  and  it  was  conjectured  that  these  differences  might 
be  caused  by  nonrandom  behavior  in  the  underlying  pseudo  random  number 
generators.  Testing  of  the  pseudo  random  number  generators  was  thus 
begun.  It  is  conjectured  that  if  the  generators  cannot  be  rejected  for 
nonrandom  behavior  using  a  set  of  standard  statistical  tests  for  random¬ 
ness,  then  semantic  differences  in  the  implementation,  instantiation, 
and/or  interpretation  of  these  two  versions  of  GPSS  are  likely.  Additional 
tests  for  these  semantic  differences  are  reported  elsewhere. 

The  following  sections  provide  a  discussion  of  the  testing  of 
the  generators,  and  the  results  of  those  tests. 

A.;?  TESTS  SELECTED 

A. 2.1  Introduction  to  Randomness  Tests. 


Three  standard  tests  of  randomness  were  chosen  in  this  study, 
namely:  (1)  the  runs  above  and  below  the  median  test,  (2)  the  maximum  of 
five  test  and  (3)  the  runs  up  and  down  test.  Each  of  these  tests  attempts 
to  determine  if  a  generated  sequence  of  numbers  is  sufficiently  random 
by  detecting  either  cyclical  patterns  or  otherwise  nonrandom  behavior. 

All  of  the  chosen  tests  are  empirical  in  that  a  computer  manipulates 
groups  of  numbers  from  the  sequence  and  computes  certain  statistics 
which  are  compared  with  standard  statistical  tablesU25l.  While  it  is 
recognized  that  there  are  a  great  many  randomness  tests,  these  particular 
tests  were  chosen  both  because  of  their  reputed  reliability  and  the  ease 
with  v/hich  their  algorithms  could  be  adapted  to  a  computer  programC25]. 
Also,  runs  tests  are  perhaps  the  only  statistical  tests  which  focus  on 
the  order  in  sequenceC26l. 


A. 2. 2  Runs  Above  and  Below  the  Median  Test 


The  first  test  chosen  was  the  runs  above  and  below  the  median 
test.  In  this  test  a  run  is  defined  as  a  series  of  either  numbers  (or 
in  the  nonparametric  approach,  ranks)  within  the  sequence  having  values 
strictly  above  or  strictly  below  the  value  of  the  median  observation. 
The  nonparametric  test  method  merely  requires  an  ordered  set  of  ranks, 
that  is,  the  relative  positions  of  the  values  of  the  observations 
within  the  sequence.  Order  is  important  because  this  test  is  based  on 
runs. 


A  test  statistic,  i.e.,  a  random  variable  whose  values  are  deter¬ 
mined  by  sample  data  [271,  can  be  calculated  based  on  the  total  number 
of  runs  in  a  sequence.  This  statistic  may  reveal  nonrandom  behavior  in 
that  either  too  few  runs  or  too  many  runs  would  likely  be  the  result  of 
a  trendy  or  cyclical  pattern.  The  sampling  distribution  of  the  number 
of  runs  can  be  approximated  by  a  normal  distribution;  therefore,  a 
normal  test  is  applied  to  the  actual  number  of  runs  in  the  sequence  [271. 

The  test  statistic  Z  is  defined  as  follows: 


Z  = 


u  -  E(u) 


» 


where 


u  =  number  of  runs  in  the  sequence, 
Zn.n, 


var(u) 


2n^n2  (Zn^ng  -  -  n^) 

(n.|  +  (n^  +  -  1)  ’ 


n]=  number  of  observations  above  the  median,  and 

n^  =  number  of  observations  below  the  median. 

The  test  statistic  Z  is  then  compared  to  critical  values  obtained  for 
a  tv/o-tailed  normal  test  from  which  the  critical  region  (the  region 
where  the  hypothesis  of  randomness  must  be  rejected)  is  determined. 

A  two-tailed  normal  test  assumes  a  normal  distribution  about  some 
mean,  and  then  a  critical  region  is  obtained  for  both  the  upper  and 
lower  tails  of  the  distribution.  If  Z  falls  within  the  critical  re¬ 
gion,  then  the  sequence  is  suspect  and  the  generator  for  the  sequence 
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is  dismissed  as  being  nonrandom.  A  negative  value  of  Z  falling  in  the 
rejection  region  implies  that  there  are  not  enough  runs  in  the  sequence; 
on  the  other  hand,  a  positive  value  of  Z  falling  in  the  rejection  region 
is  indicative  of  too  many  runs  and  possibly  a  repetitious  pattern  [27], 

A. 2. 3  Maximum  of  Five  Test. 


The  second  test  chosen  was  the  maximum  of  five  test.  Knuth  [251 
points  out  that  the  use  of  this  test  for  a  moderately  sized  sequence 
wf^l  tend  to  detect  both  local  and  global  nonrandom  behavior.  Local 
nonrandom  behavior  could  likely  be  the  result  of  clustering  of  observa¬ 
tions  around  a  single  value  while  global  nonrandom  behavior  might  be  due 
to  the  multiplier  for  the  generator  not  being  large  enough  (e.g.,  see 
Section  A. 4). 

This  test  consists  of  obtaining  observations  Hgj,  Hsj+i . 

'J5j+4  for  j  =  0,  ...,  m  -  1  where  m  is  the  integer  quotient  or  n 
divided  by  5,  n  being  the  total  number  of  observations;  let 
be  the  maximum  of  each  of  these  sequences  of  five  numbers.  The 
Kolmogorov- Smirnov  (KS)  test  method  for  measuring  the  amount  of 
deviation  between  an  assumed  distribution  function  and  the  empirical 
or  actual  distribution  function  is  used  here.  The  KS  test  is 
applied  to  the  sequence  Vg,  ....  Vq^i,  which  is  assumed  to  have  the 
cumulative  distribution  function  f(x)  =  x^  (q<x<l).  It  can  be  shown 
that  the  distribution  function  for  the  Vj's  is  indeed  F(x)  r25T.  The 
Kolmogorov-Smirnov  test  statistics  K+m  and  K-m  are  then  compared  to 
standard  statistical  tables  to  determine  whether  the  values  lie  within 
the  critical  regions  for  given  confidence  levels,  where  K+m  is  the 
greatest  amount  of  deviation  vfhen  the  actual  distribution  function  is 
greater  than  F(x),  and  K“m  is  the  greatest  amount  of  deviation  when 
the  actual  distribution  function  is  less  than  F(x).  If  the  values  of 
K*m  or  K"m  are  in  the  critical  regions,  then  the  hypothesis  that  the 
sequence  is  random  must  be  rejected. 

A. 2. 4  Runs  Up  and  Down  Test. 

The  last  of  the  three  tests  selected  was  the  runs  up  and 
down  test.  This  test  is  examined  in  detail  by  both  Knuth  [25]  and 
Fishman  [287.  The  associated  test  statistic  is  calculated  based  on 
the  number  of  runs  up  and  the  number  of  runs  dov<n.  Here,  a  run  is 
defined  as  a  series  of  observations  such  that  <...<  X^+^ 

for  runs  up,  or  conversely,  Xj  >  Xj+j  >...>  Xj+j  for  runs  down, 
for- r,s>o.  The  test  statistic  is  given  by 


P  P 

"  ’  j?i  Sj  tKj-  E{Rj)] 


where  »  number  of  runs  of  length  1, 

Rj  ■  number  of  runs  of  length  .1, 

E(Rfj  =  expected  number  of  runs  of  length  1  (see  Table  A.l), 

E{R.1)  *  expected  number  of  runs  of  length  J  (see  Table  A.l), 

Cjj  =  element  In  row  1  and  column  j  of  the  Inverse  of  the 
covariance  matrix  of  Rj,  ....  Rp  (see  Table  A. 2), 
p  =  length  of  longest  run. 

TABLE  A.l  (from  Fishman  [2ftl) 

E(R,)  .  1  -.4  . 

’  (1  +  3)f  (1  +  3)1 


=0.4167n  + 

0.0933 

1=1 

=n.l933n  - 

0.2333 

1=2 

=n.n52fin  - 

0.1306 

1=3 

=n.0115n  - 

0.0413 

1=4 

=0.0n20n  - 

0.0095 

1=5 

=0.00n3n 

0.0017 

1=6 

*3.9xin-5n  . 

0.0003 

1=7 

TABLE  A.2 

4529.4 

9044.9  13569 

1P091 

22615 

27992 

9044.9 

19097  27139 

36197 

45234 

55799 

13568 

27139  40721 

54291 

67952 

93695 

lanoi 

36197  54291 

72414 

90470 

111590 

22615 

45234  67952 

90470 

113262 

139476 

27992 

55799  83695 

111590 

139476 

172860 

The  test  statistic  R  Is  known  to  have  an  asymptotically  chi-square  (11 
tributlon  with  p  degrees  of  freedom  r2B!l.  Fishman  proposes  an  analo¬ 
gous  form  using  a  six  degree  of  free<jom  chi -square  distribution  for 
either  of  the  cases  v/here  p  =  *5  or  p  »  7.  This  form  combines  Ry  with 
Rp  and  ElRy)  with  E(Rp).  When  p  Is  equal  to  five,  Rp  Is  set  to  zero 
so  that  the  computer  program  used  for  the  testing  need  not  he  altered 
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A. 3  TESTING 


A. 3.1  Judgment  Criteria. 

For  the  analysis  of  the  “goodness”  of  a  pseudo  random  number 
generator,  the  criteria  given  by  Knuth  [?5]  were  used.  The  criteria 
specify  that  for  the  range  of  a  given  statistic  S,  a  generator  Is  classi¬ 
fied  as  rejected  if  the  value  computed  for  a  sample,  S*,  lies  In  the 
outermost  two  percent  of  the  known  distribution  function  of  S  (one  percent 
on  each  end).  Likewise,  it  Is  classified  as  "suspect"  If  S*  lies  In  the 
next  Innermost  eight  percent  and  "almost  suspect"  If  It  lies  In  the  next 
Innermost  ten  percent.  The  following  table  summarizes  these  criteria. 

TABLE  A.3.  ACCEPTANCE  INOICATORS  VERSUS  TEST  STATISTICS 


S*  In  Range  of  S 
0-1  percent,  99-100  percent 
1-5  percent,  95-99  percent 
5-10  percent,  90-95  percent 


Indication 

Reject 

Suspect 

Almost  Suspect 


Translating  this  table  to  the  particular  tests  being  used  gives  crit¬ 
ical  regions  as  shotvn  In  Table  A,4. 

One  additional  consideration  should  be  examined  concerning 
the  use  of  multiple  tests.  For  a  rejection  region  of  size  alpha 
using  N  tests,  the  probability  of  rejecting  a  generator  even  though 
the  hypothesis  of  randomness  Is  true  Is  given  by  1  -  (1-alpha)^.  Here 
alpha  *  0.02  and  N  =  3,  so  the  probability  of  rejecting  a  generator  that 
Is  actually  random  enough  Is  1  -  (1-0.02)^  =  0.06;  therefore,  the  cri¬ 
teria  of  rejection  used  In  this  study  lead  to  a  94  percent  confidence 
level. 


A. 3. 2  Test  Procedures. 

The  ten  UNIVAC  RPSS  1100  pseudo  random  number  generators  were 
tested  along  with  that  of  I5W  GPSS/360.  GPSS/360  actually  has  eight 
generators  available,  but  when  they  are  used  in  unmodified  form,  each 
returns  an  identical  sequence  of  random  numbers  [1,  p.lAAl.  The  UNIVAC 
generators  were  tested  by  using  the  GPSS  1100  random  number  generation 
algorithm  to  produce  a  sequence  of  numbers.  Using  the  algorithm.  In¬ 
stead  of  merely  copying  a  sequence  of  numbers  from  a  GPSS  1100  program, 
saved  considerable  time.  It  should  be  noted  that  the  sequence  generated 
by  this  approach  was  checked  against  the  output  of  actual  random  numbers 
from  a  GPSS  1100  program  to  Insure  exact  replication  of  the  sequences. 
Unfortunately,  this  approach  could  not  be  easily  applied  to  the  IBM 
generator.  This  prompted  the  writing  of  a  short  program  In  GPSS/360  In 
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TABLE  A.4.  ACCEPTANCE  INDICATORS  VERSUS  TEST  STATISTIC  CRITICAL  REGIONS 


1.  Runs  above  and  below  the  median 
|Z|  >  2.33 

2.33  >  |Z|  >  1.65 
1.65  >  |Z|  >.  1.28 

2.  Maximum  of  five 

K600  <  0.0648 


Reject 

Suspect 

Almost  Suspect 

Reject 


K600  I  1.5092 
0.648  <  K600  <  .1544 

1.5092  >  K600  >_  1.2170 
K2n0  <  .0603 

K200  1.5033 

0.0603  <  K200  <  .1502 

1.5033  >  K200  >  1.2119 

3.  Runs  up  and  down 
R  <  .872 

R  >_  16.81 

.872  <  R  <  1.64 
16.81  >  R  12.59 

1.64  <  R  <  2.20 

12.59  >  R  >  10.65 


Suspect 

Reject 

Suspect 

Reject 

Suspect 

Almost  Suspect 
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order  to  provide  a  listing  of  the  IBM  sequence,  which  was  then  read  into 
the  testing  program.  Only  the  results  of  tests  for  sequences  of  length 
1000  to  3000  are  discussed  in  detail  because  the  simulation  models  of 
concern  in  this  study  call  on  any  given  generator  approximately  that  many 
times  in  any  run.  Tests  on  sequences  of  length  greater  than  5000  are  of 
little  interest  at  this  point,  but  some  results  of  tests  on  these  longer 
sequences  are  given  in  Table  A, 6,  A  discussion  of  the  random  number  gener¬ 
ation  techniques  is  given  in  the  next  section, 

A, 4  GPSS  PSEllOO  RANDOM  NUMBER  GENERATION  SCHEMES 


A,4,l  IBM, 

According  to  the  IBM  GPSS/360  User’s  Manual  [2,  pp,  36-371,  the 
random  number  generation  algorithm  is  as  follows: 

1,  The  appropriate  word  of  the  index  array  points  to  one  of  the 
eight  numbers  in  the  base  number  array.  Since  the  index  array  words  are 
initially  zero,  the  first  base  number  used  will  be  the  seed, 

2,  The  appropriate  number  in  the  multiplier  array  is  multi¬ 
plied  by  the  base  number  chosen  in  step  1, 

3,  The  low-order  31  bits  of  this  product  are  stored  in  the 
appropriate  word  of  the  multiplier  array,  to  be  used  the  next  time  a  random 
number  is  called  for, 

4,  Three  bits  of  the  high-order  16  bits  of  the  product  pro¬ 
duced  in  step  2  are  stored  in  the  appropriate  word  of  the  index  array,  for 
future  use.  This  number  (0-7)  points  to  one  Of  eight  words  of  the  base 
number  array  to  be  used  the  next  time  a  random  number  is  called  for, 

5,  (a)  If  the  random  number  required  is  a  fraction,  the  middle 
32  bits  of  the  product  produced  in  step  2  are  divided  by  lO^,  and  the 
remainder  becomes  the  six-digit  fractional  random  number, 

(b)  If  the  random  number  required  is  an  integer,  th;?  middle 
32  bits  of  the  product  produced  in  step  2  are  divided  by  l63,  and  the 
remainder  becomes  the  three-digit  random  number, 

A,4,2  UN  IV AC, 

The  UNIVAC  random  number  generation  algorithm  [3,  pp,3,30,  3,32] 
is  a  sinple  one.  It  uses  a  linear  congruential  or  mixed  linear  congruential 
generator,  as  the  case  may  be.  It  takes  the  form 

X  j  ®  s 

Xn+1  =  (mXp+I)  mod  2^5 

where  S*  seed,  m»  multiplier,  and  I*  incrercnt.  When  a  fractional  number 
is  needed,  the  integer  is  divided  by  When  an  Integer  value 
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from  n  to  999  is  required,  the  fractional  number  is  multiplied  by  103 
and  truncated. 

A. 4. 3  Independent  Streams  of  Random  Numbers. 

The  UNIVAC  pseudo  random  number  generator  uses  ten  different 
combinations  of  multipliers,  increfnents,  and  seeds  to  produce  its  ten 
random  number  sequences.  The  IBM  has  one  generator,  replicated  eight 
times. 

A. 5  RESULTS  OF  TESTS 

Using  the  established  critical  regions,  it  can  be  seen  that 
most  of  the  generators  fared  well.  (See  Table  A. 5.)  It  appears  that 
UNIVAC  generator  nine  may  have  a  few  problems  associated  with  its  use; 
the  values  of  the  runs  up  and  dovm  test  statistics  for  sequence  sizes  of 
both  1000  and  3000  lie  in  the  rejection  region.  Also,  the  value  of  the 
maximum  of  five  test  statistic  K“600  places  more  suspicion  on  the  se¬ 
quence  produced  by  this  generator.  These  facts  suggest  that  generator 
nine  should  not  be  used,  at  least  in  short  simulation  models,  because 
the  number  sequence  produced  by  it  does  not  exhibit  sufficient  randomness. 

The  only  other  generators  with  test  statistic  values  in  the 
rejection  region  are  the  UNIVAC  generators  one  and  two.  The  runs  up  and 
down  test  statistic  for  a  sequence  length  of  1000  is  far  too  large  for 
each  of  the  generators.  It  is  interesting  to  note  that  generator  one  is 
used  as  the  resident  generator  in  the  OPSS  1100  language.  This  means 
that  on  occasions  when  the  TIME  and  00  TO  fields  require  a  random  number, 
they  call  on  generator  one.  (It  should  also  he  noted  that  the  simulation 
models  studied  did  not  include  these  types  of  TIME  and  GO  TO  fields.) 
Generator  two,  which  should  also  be  rejected  for  a  sequence  size  of  1000 
according  to  Knuth's  criteria,  was  employed  in  all  four  of  the  UNIVAC 
simulation  models  studied.  For  each  message  introduced  into  the  model, 
the  generator  was  called  on  twice,  once  to  generate  Poisson  interarrivals, 
and  once  to  create  exponentially  distributed  message  lengths.  Since  a 
minimum  of  1000  messages  were  included  in  each  run,  the  second  generator 
was  called  on  at  least  2000  times,  probably  closer  to  3000  times  when 
"warmup"  and  queued  messages  are  counted.  Therefore,  the  nonrandom 
behavior  of  generator  two  for  a  sequence  size  of  1000  does  not  appear  to 
he  a  possible  cause  for  the  discrepancy  between  the  UNIVAC  and  IBM  results. 

The  only  other  generator  that  is  reasonably  suspicious  is  the 
third  UNIVAC  generator.  Three  of  the  four  maximum  of  five  test  statistics 
for  sequences  from  this  generator  lie  in  the  "suspicion"  range.  Inciden¬ 
tally,  this  is  the  generator  used  in  the  uniform  distribution  function  in 
the  UNIVAC  models  used  to  determine  the  routing  of  the  messages. 

Since  only  two  random  number  generators  are  required  for  the 
UNIVAC  simulation  models  in  addition  to  generator  one,  it  would  seem 
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advantageous  to  select  generators  that  cast  the  least  doubt  on  the  re> 
suits.  This  usage  of  the  "best"  generators  would  lead  to  a  more  meaning¬ 
ful  comparison  between  UN IV AC  and  IPM  data. 

There  appears  to  be  no  need  to  tamper  with  the  IBM  generator 
as  It  comes  through  the  tests  very  well.  Rut,  If  longer  sequences  are 
accepted  for  UNIVAC,  then  IBM  sequences  of  similar  length  should  be 
tested  for  randomness. 

Examination  of  even  longer  sequences  for  the  UNIVAC  generators 
(see  Table  A.6)  shows  a  trend  for  almost  all  of  the  generators  falling 
the  runs  above  and  below  the  median  test  for  sequence  sizes  greater  than 
10,000  numbers.  The  maximum  of  five  test  and  the  runs  up  and  down  test 
reject  generators  seven  and  six,  respectively,  for  sequences  of  8000 
numbers  and  up.  From  these  results.  It  can  he  seen  that  there  are  partic¬ 
ular  generators  that  should  be  avoided  for  certain  sequence  sizes. 

A.6  SUMMARY  AND  CONCLUSION 

The  randomness  tests  performed  Indicate  that  the  UNIVAC  gen¬ 
erators  are  primarily  suited  for  models  requiring  numerical  sequences 
of  length  from  3000  to  somewhere  around  8000.  The  IBM  generator  cannot 
be  rejected  at  the  94  percent  confidence  level  for  sequences  of  length 
1000  or  3000,  but  a  study  of  Its  characteristics  for  longer  sequences 
should  be  performed.  From  this  study.  It  appears  the  generators 
used  In  the  simulation  models  are  In  fact  random  enough  and  do  not  cause 
the  principal  differences  between  UNIVAC  and  IBM  simulation  results. 
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TABLE  A. 5  SUMMARY  OF  RAHOOM  NUMBER  TESTS 


HENERATOR 

SEQUENCE 

SHE 

MEAN 

MEOIAM 

K+(N/5) 

K-(N/5) 

R(M) 

UNIVAC  1 

1000 

SOS 

513 

0.13 

0.7739 

0.3101 

31.74*** 

3000 

S02 

509 

-0.22 

0.8147 

0.4625 

7.22 

UHIVAC  2 

1000 

4«4 

475 

0.63 

0.7957 

0.1212** 

33.15*** 

3000 

406 

492 

0.07 

0.0345 

0.1802 

10.52 

UMIVAC  3 

1000 

407 

400 

-1.27 

0.9341  ^ 

0.0891** 

2.94 

3000 

40S 

402 

-1.02 

1.3372** 

0.0727** 

7.36 

UMIVAC  4 

1000 

40fl 

40Q 

-0.25 

0.8428 

0.1650 

15.81** 

3000 

401 

4P7 

0.04 

0.9271 

0.6819 

1.17** 

UNIVAC  S 

1000 

Aop 

4BB 

-1.71** 

0.8114 

0.2349 

4.41 

3000 

sio 

511 

-1.04** 

0.3244 

1.1122 

5.24 

UNIVAC  fi 

1000 

506 

509 

1.39* 

0.2052 

0.8640 

6.68 

3000 

402 

491 

1.50* 

0.7937 

0.2979 

2.R2 

UNIVAC  7 

1000 

4ft4 

469 

0.76 

1.1430 

0.3604 

5.00 

3000 

490 

406 

1.17 

0.8742 

0.6500 

5.48 

UNIVAC  B 

1000 

409 

514 

-0.70 

0.8881 

0.5400 

12.71** 

3000 

496 

494 

0.07 

0.5209 

1.0026 

2.99 

UNIVAC  0 

1000 

516 

517 

0.00 

0.5800 

0.9238 

34.12*** 

3000 

501 

50fl 

0.84 

0.3287 

1.3742** 

19.69*** 

UNIVAC  10 

1000 

502 

505 

-1.90** 

0.5561 

1.0451 

4.33 

3000 

405 

498 

-1.20 

0.9188 

0.6737 

2.19* 

IBM 

1000 

497 

484 

-1.45* 

0.9675 

0.3310 

11.07* 

3000 

493 

490 

0.44 

1.0176 

0.1103** 

7.55 

*  Alnost  suspect 

**  Suspect 

***  Reject 


TABLE  A.6  SIW^ARY  Of  BANOW  NUMBER  TESTS 


SF^fE!’ATOR 

NUMBER  OF 
RENERATEO 
R.N.'S 

BASIC  SERIES 
STATISTICS 

RUNS  ABOVE 
Afin 

BELOV)  TUE 
MEDIAN  TEST 

MAX  IWJM  OF  5  TEST 

RUNS  UP  AND 
DOWN  TEST 

NUMBER 

M 

MEAN  MFniAN 

Zn 

K+(N/5)  K-(N/5) 

R(N) 

APPENDIX  B 

GPSS/360  PR06RA«1  LISTINGS 
FOR  RING  NETWORK  SIMULATIONS 
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Next  page  Is  blank. 


Fo*'  NEWHALL/IBM  GPSS  Program  Listing  see  Reames  [6]»  pp.  191-194. 

The  following  blocks  were  Inserted  at  the  top  of  the  program  shown  In 
Reames  [6]  In  order  to  successfully  execute  the  GPSS/360  program  on 
the  APG  IBM  360/65  computer  system: 


REALLOCATE  XAC,1200.BLO,100,FAC.100,ST0.100,QUE.100.L06,100 
REALLOCATE  TAB . 50 . FUN . 1 0 . VAR . 20 . FSV , 1 00 . HSV . 50 , CHA . 1 00 
REALLOCATE  BVR .1 O.FMS .1 O.HMS .1 0 .MAC. 5 .COM.90000 

SIMULATE 


For  the  PIERCE/IBM  program  listing  see  Reames  [6],  pp.  187-190. 

The  change  to  this  program  starts  at  the  bottom  of  page  189  In  [6] 
and  Is  as  follows: 


*  LAST  PACKET  OF  A  MESSAGE  HAS  BEEN  RECEIVED.  RECORD  TOTAL 

*  MESSAGE  TRANSMISSION  TIME. 

* 

LASTP  TABULATE  TMGTM  RECORD  TOTAL  MESSAGE  TRANSIT  TIME 

* 

*  CHECK  IF  LAST  TERMINATION  THEN  SAVE  RELATIVE  CLOCK 

* 

SAVEVALUE  3+,Kl 
TEST  E  X3,X4,PATW 
SAVEVALUE  2+.C1 

* 

*  SAVES  VALUE  OF  RELATIVE  CLOCK  FOR  ABSOLUTE  CLOCK 

* 

PATW  TERMINATE  1 

* 


*  TABLES  AND  QTABLES  — 


0 


For  DLCNNE/IBM  Program  Listing  see  Reames  [6],  pp.  178-186. 

DLCNNE  Is  Identical  to  DLCN  except  that  the  following  blocks  In  DLCN  at 
the  top  of  page  181  In  [6],  which  now  read: 


RECVR 

LOGIC  S 

*1 

TRANSFER 

.010,*+4.*+l 

TRANSFER 

.010.*+3,*+l 

ASSIGN 

5,K3 

TRANSFER 

.RECVD 

LOOP 

6,RECVR+1 

RECVD 

ADVANCE 

V$AMSG 

have  been  changed 

In  DLCNNE  to  re 

RECVR 

LOGIC  S 

..  ^ 

TRANSFER 

.*+3 

ASSIGN 

5,K3 

TRANSFER 

.RECVD 

ft 

LOOP 

6,RECVR+1 

RECVD 

ADVANCE 

V$AMSG 

GET  CONTROL  OF  RECEIVER 
PERFORM  MSG  ERROR  CHECKING. 
ASSUMING  1  ERROR  PER  10,000  CHARS. 
IF  ERROR.  SET  ACK  MSG  RESPONSE 
&  GO  SEND  ACK  MSG 
CHECK  EACH  CHAR.  OF  MSG  FOR  ERROR 

ALLOW  TIME  TO  RECEIVE  DATA 


GET  CONTROL  OF  RECEIVER 
SKIP  POSSIBILITY  OF  ERRORS  IN  CHARS. 
IF  ERROR.  SET  ACK  MSG  RESPONSE 
&  GO  SEND  ACK  MSG. 

CHECK  EACH  CHAR.  OF  MSG  FOR  ERROR 

ALLOW  TIME  TO  RECEIVE  DATA 


This  change  disables  retransmissions  due  to  received  character  errors; 
hence,  the  name  DLCN/"Np  Errors"  or  simply  DLCNNE.. 
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APPENDIX  0 


GLOSSARY 


AMSAA 

APG 

ARRADCOH  - 
ASAS  FSD  - 

c3a 

COC 

CPU 

CSD 

DA 

DLCN 

DLCNNE 

GPSS 

IBM 

OPTADS 

PM 

SACDIN 

SIGMA 

SIMSCRIPT  - 


TOS  CASE 
UNIVAC 


US  Am^  Materiel  Systems  Analysis  Activity 

Aberdeen  Proving  Ground,  Maryland 

US  Army  Armament  Research  and  Development  Command 

All  Source  Analysis  System  Full  Scale  Development 

Command,  Control,  and  Communications  Analysis 

Trademark  and  abbreviation  for  the  Control  Data  Corporation 

Central  Processing  Unit  of  computer  systems 

Combat  Support  Division 

Department  of  the  Army 

Distributed  Loop  Computer  Network,  The  Ohio  State  University. 

Modified  simulation  model  for  DLCN  with  no  errors  In  character 
transmi sslon 

Either  of  two  simulation  language  dialects  called  "General 
Purpose  Simulation  System"  by  IBM  and  called  “General  Purpose 
Systems  Simulator"  by  UNIVAC 

Trademark  and  abbreviation  for  International  Business  Machines 
Corporation 

Operations  Tactical  Data  ^sterns 

Program  or  Project  Manager 

Stragetic  Air  Command  Digital  Network 

Name  of  force  level  maneuver  control  system 

Generic  name  of  a  computer  programming  language  developed  at 
the  RAND  Corporation  for  discrete  event  simulation  with  a 
version  marketed  under  the  trademark  SIMSCRIPT  1 1.5  by  Consoli¬ 
dated  Analysis  Centers,  Inc. 

Tactical  Operations  Systems  for  Corps  and  Subordinate  Echelons 

Trademark  and  name  of  the  Sperry  UNIVAC  Division  of  the  Sperry 
Rand  Corporation 
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